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Abstract 


Every  littoral  or  regional  ocean  has  both  coastal  circulation  components  and  influences  from 
offshore  regions.  Over  the  past  decade,  the  feature-oriented  regional  modeling  approach  has 
evolved  to  become  operational  in  some  of  these  regions  (western  North  Atlantic,  tropical  North 
Atlantic).  This  methodology  is  model-independent  and  can  be  applied  with  or  without  satellite 
and/or  in  situ  observations,  especially  in  coastal  regions.  In  applying  this  methodology  for  the 
subtropical  North  Pacific,  we  first  present  an  overview  of  the  synoptic  nature  of  the  different 
features  that  comprise  the  multiscale  and  complex  circulation  in  the  California  Current  system 
(CCS).  This  description  is  a  prerequisite  to  feature  modeling  and  is  based  on  our  present 
understanding  of  the  features  and  their  dominant  space  and  time  scales  of  variability.  The  region 
of  interest  is  from  36°N  to  48°N  -  essentially  most  of  the  eastern  boundary  current  region  of  the 
subtropical  North  Pacific.  We  present  a  synergistic  configuration  of  synoptic  features  interacting 
with  one  another  on  multiple  and  sometimes  overlapping  space  and  time  scales  as  a  meander- 
eddy-upwelling  system.  The  second  step  is  to  define  the  feature-oriented  regional  modeling 
system  (FORMS)  for  this  system.  The  major  multiscale  synoptic  circulation  features  include  the 
mean  flow  and  southeastward  meandering  jet(s)  of  the  California  Current,  the  poleward  flowing 
California  Undercurrent,  and  six  upwelling  regions  along  the  western  coast  of  the  U.S.  We  first 
identify  the  typical  synoptic  width,  location,  vertical  extent,  and  core  characteristics  of  these 
features  and  their  dominant  spatial  and  temporal  scales  of  variability  from  past  observational, 
theoretical  and  modeling  studies.  These  synoptic  features  are  then  melded  with  available 
climatology  of  the  region  for  initialization  and  assimilation. 

Finally,  dynamical  simulations  are  run  as  nowcasts  and  short-term  (4-6  weeks)  forecasts 

using  these  feature  models  as  initial  fields  and  the  Princeton  Ocean  Model  (POM)  for  dynamics. 
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The  set  of  simulations  over  a  40-day  period  illustrate  the  applicability  of  FORMS  to  a  transient 
eastern  boundary  current  region  such  as  the  CCS.  The  success  of  the  dynamical  simulations  were 
realized  in  (i)  maintenance  of  the  low-salinity  pool  in  the  core  of  the  CC;  (ii)  representation  of 
eddy  activity  inshore  of  the  coastal  transition  zone;  (iii)  realistic  eddy  kinetic  energy  evolution; 
(iv)  subsurface  (intermediate  depth)  mesoscale  feature  evolution;  and  (v)  deep  poleward  flow 
evolution.  (Keywords:  California  Current  System;  Feature  Models;  FORMS;  Upwelling) 

1.  Introduction 

The  California  Current  system  (CCS)  is  at  the  eastern  boundary  of  the  wind-driven 
subtropical  gyre  in  the  North  Pacific.  The  region  of  interest  for  this  multiscale  (in  space  and 
time)  overview  is  from  36°N  to  48°N  and  from  145°W  to  the  North  American  west  coast.  In 
short,  large-scale  gyres  and  mean  southeastward  flow  are  prevalent  in  climatological  averages, 
while  mesoscale  eddies  and  submesoscale  filaments  are  abundant  in  its  synoptic  state  and  are 
transient  in  nature.  The  regional  ecosystem  is  well  known  for  its  productivity  and  the  seasonal 
and  inter-annual  variability  resulting  from  large-scale  wind-driven  upwelling  around  its  multiple 
capes  and  their  interactions  with  the  irregular  occurrences  of  the  El  Nino  -  La  Nina  phenomenon. 
The  typical  summer  circulation  is  drastically  different  from  that  in  the  winter,  and  a  synoptic 
description  at  any  one  time  has  been  a  challenge. 

This  region  has  been  studied  systematically  from  the  early  1950s  by  the  CALCOFI 
(California  Cooperative  Oceanic  Fisheries  Investigation;  Wyllie,  1966)  monitoring  program. 
However,  as  in  many  other  regions,  the  richness  of  mesoscale  variability  was  recognized  first  in 
satellite  SST  images  for  this  region,  and  then  in  color  images  now  available  from  MODIS  and 
other  satellite  sensors.  Figure  1  shows  the  difference  in  features,  their  appearance,  and  prevalent 
scales  as  available  from  an  SST  image  from  those  of  a  typical  climatology.  Clearly,  the 
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averaging  of  data  at  larger  scales  masks  the  presence  of  many  of  the  important  mesoscale  and 
submesoscale  features  of  circulation. 

This  region  of  interest  —  i.e.,  the  western  coast  of  the  U.S.  —  includes  both  a  larger-scale 
dynamic  offshore  region  and  a  finer-scale  dynamic  coastal  region.  The  offshore  region  is 
primarily  dominated  by  the  large-scale  California  Current,  the  California  Undercurrent  and  parts 
of  subtropical  and  sub-polar  gyre  circulations  in  the  eastern  Pacific.  The  coastal  region  includes 
features  such  as  upwelling  fronts,  cold  pools  inshore  of  these  fronts,  filaments,  squirts, 
mushroom-head  vortices,  mesoscale  and  sub-mesoscale  eddies,  and  meanders. 

One  of  the  approaches  for  regional  modeling  of  such  regions  in  the  world  oceans  is  the 
use  of  “knowledge-based  feature  models.”  This  approach  is  distinctly  different  from  the  basin- 
scale  approach  that  requires  models  to  develop  the  inertia  fields  in  a  so-called  "spin-up"  period, 
which  could  be  anywhere  from  one  to  ten  years.  Our  goal  is  to  develop  a  system  for  short-term 
synoptic  forecasting  for  use  where  we  cannot  afford  to  run  models  for  a  long  period  of  time  for 
spin-up  or  if  model  biases  could  build  up  when  numerical  codes  are  run  for  long  periods.  This 
feature-oriented  approach  has  been  used  for  regional  simulations  and  operational  forecasting  for 
the  past  two  decades.  It  originated  and  was  very  successful  in  weather  prediction  when  sufficient 
data  were  not  available  for  a  complete  initialization  of  features  (e.g.  Bennett,  1992;  Bennett, 
2002).  In  the  ocean,  specific  examples  include  the  studies  by  Robinson  et  al.  (1989),  Hurlburt  et 
al.  (1990,  1996),  Fox  et  al.  (1992),  Glenn  and  Robinson  (1995),  Cummings  et  al.  (1997), 
Gangopadhyay  et  al.  (1997),  Lermusiaux  (1999a),  and  Robinson  and  Glenn  (1999)  that  have 
applied  the  feature  modeling  technique  for  use  in  nowcasting,  forecasting  and  assimilation  of 
various  in  situ  (XBT,  CTD)  and  satellite  (SST,  SSH  [GEOSAT,  TOPEX/Poseidon]  and  SSC 
[SeaWiES])  observations  in  the  western  North  Atlantic. 
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Gangopadhyay  and  Robinson  (2002)  have  generalized  the  feature-oriented  approach  for 
strategic  application  to  any  regional  ocean.  A  feature-oriented  regional  modeling  system 
(FORMS)  for  the  Gulf  of  Maine  and  Georges  Bank  region  has  been  developed  for  real-time 
applications  for  medium-range  (7-10  days)  and  mesoscale  to  sub-mesoscale  forecasting 
(Gangopadhyay  et  al.,  2003;  Brown  et  al.  2007a,b).  This  feature-oriented  methodology  is  also 
model-independent  and  can  be  applied  in  lieu  of  either  satellite  or  in  situ  observations,  or  both, 
especially  in  coastal  regions.  It  has  been  applied  in  many  real-time  simulations,  including  the 
recent  lOOS  efforts  (Schofield  et  al,  2010). 

Feature-oriented  methodology  requires  developing  a  synoptic  circulation  template  for  the 
region,  called  the  “basis  template.”  This  template  is  designed  from  a  feature-based  synthesis  of 
the  regional  circulation  patterns.  From  this  “basis  template,”  a  map  of  strategic  sampling 
locations  for  placing  feature  model  profdes  is  produced.  These  profdes  provide  “synthetic 
synoptic  expressions”  for  fronts,  jets,  eddies,  gyres,  and  other  circulation  structures  and  water 
masses  at  the  initialization  or  updating  phases.  Features  are  represented  by  both  analytical 
structures  and  by  synoptic  sections  and  profiles.  The  feature  model  profiles  for  a  typical  feature 
are  developed  by  analyzing  past  synoptic  high-resolution  observations  for  that  feature. 
Gangopadhyay  and  Robinson  (2002)  described  the  generalized  mathematical  forms  for  some  of 
the  typical  feature  models  developed  so  far  (see  their  Appendix  A).  Gangopadhyay  et  al.  (2003) 
described  a  number  of  applications  for  coastal  regions  such  as  a  coastal  current,  a  tidal  mixing 
front,  eddies  and  gyres. 

For  any  oceanic  region  such  as  the  California  Current  system,  such  structures  and 
parameterized  forms  can  be  implemented  to  characterize  the  relevant  circulation  entities.  Such 
implementation  requires  careful  and  detailed  scientific  analyses  to  identify  the  spatial  and 
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temporal  scales  and  variability  that  define  and  distinguish  these  features  from  one  another  while 
preserving  the  particular  characteristics  and  dynamic  balances  of  each  feature.  After  major 
features  are  identified  and  individually  represented,  they  are  used  in  the  initialization  of  a  basic 
dynamical  model  (e.g.,  POM).  Dynamical  adjustment  accomplishes  two  important  tasks:  i)  a 
consistent  dynamical  interaction  of  the  features,  and  ii)  the  generation  of  smaller-scale  features, 
such  as  squirts  and  sub-mesoscale  eddies. 

In  this  study,  we  build  on  the  experience  and  knowledge  from  the  above  studies  to 
develop  FORMS  for  the  California  Current  system.  The  application  of  FORMS  for  nowcasting 
and  forecasting  is  of  interest  to  the  US  Navy,  CeNCOOS,  SCCOOS  and  similar  entities.  In 
Section  2,  we  first  provide  an  overview  of  the  synopticity  of  the  CCS  circulation.  In  Section  3, 
we  then  present  a  synopsis  of  the  CCS  features  and  their  variability,  which  is  the  prerequisite 
first  step  of  a  feature-oriented  regional  modeling  perspective.  The  individual  feature  models  are 
described  in  Section  4.  The  FORMS  protocol  for  initialization  for  this  region  is  then  illustrated  in 
Section  5.  Dynamical  simulations  using  climatology  and  selected  feature  models  are  discussed  in 
Section  6.  Section  7  summarizes  the  results,  and  also  outlines  the  future  directions  for 
maintaining  longer-term  synoptic  forecasts. 

2.  A  Brief  Overview  of  the  Synoptic  Nature  of  the  CCS 

Moores  and  Robinson  (1984)  first  described  the  instantaneous  California  Current  (CC)  as 
an  energetic  meandering  jet,  interacting  with  mesoscale  eddies  and  coastal  upwelling  regions. 
During  the  OPTOMA  (Ocean  Prediction  through  Observation,  Modeling  and  Analysis)  exercises 
in  the  summers  of  1983  and  1984,  Robinson  et  al.  (1984,  1986)  designed  the  first,  real-time 
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dynamical  forecasting  experiment  in  the  CCS  region  and  found  that  the  time-scales  of  eddy- 
meander/  eddy-eddy  interactions  are  on  the  order  of  1-2  weeks  in  this  region. 

Brink  et  al.  (1991)  defined  a  coastal  transition  zone  (CTZ)  along  the  California  coast  as  the 
demarcation  between  the  faster,  more  energetic,  shallow  and  narrow  coastal  upwelling- 
dominated  flow  system  and  the  large-scale,  broad  mean  southward  flow  in  the  deep  offshore 
region.  The  transition  zone  was  thought  to  be  an  area  where  mesoscale  eddies  and  filaments 
would  dominate.  Strub  et  al.  (1991)  proposed  a  synoptic  strong  surface  jet  as  the  CC  core 
meandering  through  a  field  of  associated  cyclonic  and  anticyclonic  eddies.  They  also  found 
evidence  of  some  subsurface  eddies  which  had  no  identifiable  surface  signatures.  Hickey  (1998) 
provided  an  end-of-the-century  overview  of  the  CCS,  its  seasonality  and  inter-annual  variability. 
Hickey  identified  the  core  of  the  current  as  a  freshwater  signature,  different  from  the  coastal 
upwelling  water  masses. 

A  compilation  of  the  physical,  chemical  and  biological  studies  of  the  CCS  in  the  past 
century  was  organized  by  Chavez  and  Collins  (2000)  in  a  special  issue  of  DEEP  SEA 
RESEARCH.  The  mean  structure  of  the  poleward  California  Undercurrent  (CUC)  was  described 
by  Collins  et  al.  (2000)  along  Pt.  Sur.  Pierce  et  al.  (2000)  showed  the  existence  and  continuation 
of  this  flow  along  the  California  coast,  while  Noble  and  Ramp  (2000)  analyzed  the  spatial  and 
temporal  variability  in  the  CUC.  The  latter  study  showed  how  and  where  the  poleward  flow  is 
surface-intensified  and  how  the  variability  is  related  to  the  surface  wind  stress.  Barth  et  al. 

(2000)  described  how  the  surface  upwelling  jets  separate  from  the  cape-like  coast  and  become  a 
conduit  for  transferring  fluxes  across  the  shelf. 

Strub  and  James  (2000)  investigated  the  seasonality  of  the  mesoscale  evolution  of  the 
CCS  based  on  altimeter-derived  surface  velocity  fields.  Eddy  statistics  such  as  eddy  kinetic 
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energy  (EKE),  wavelengths  of  meanders  and  jets  and  their  spatial  extents  were,  for  the  first  time, 
quantified.  Strub  and  James  (2000)  raised  two  very  important  questions:  one  related  to  the 
formation  and  westward  migration  of  the  seasonal  jet  and  another  related  to  the  source  of 
freshwater  within  the  CC  core.  They  hypothesized  that  the  seasonal  jet  begins  as  one  or  more 
local  upwelling-front  jets  which  move  offshore  by  Rossby  waves.  These  jets  unite  offshore  (to 
about  130W)  to  become  a  free,  open-ocean  jet  that  maintains  its  identity  as  the  CC  core  during 
spring  and  summer.  Based  on  previous  observation  and  models,  Strub  and  James  (2000)  also 
hypothesized  that  the  source  of  the  fresh  water  in  the  core  of  the  CCS  has  multiple  contributions 
from  the  North  Pacific  Current,  coastal  inputs  from  the  Columbia  River,  the  Strait  of  Juan  da 
Euca  and  the  British  Columbia  coast  (sub-arctic  origin). 

In  a  high-resolution  analysis  of  hydrographic  data  along  line  67  during  1988-2002, 
Collins  et  al.  (2003)  observed  the  California  Current  jet  at  about  100-200  km  offshore  near 
Monterey  Bay.  Their  analysis  suggests  that  the  “mean”  CC  jet  is  at  the  inshore  edge  of  the 
broader  mean  seasonal  equatorward  flow  of  the  California  Current.  The  mean  jet's  maximum 
velocity  (6-10  cm/sec)  is  much  weaker  than  that  of  the  baroclinic  jet  (50-100  cm/sec)  identified 
by  Brink  et  al.  (1991).  One  interpretation  is  that  the  average  currents  over  the  50  km  or  so  mean 
CC  jet  width  are  on  the  order  of  10  cm/sec,  but  that  much  higher  velocities  are  found  in  the 
narrow  frontal  regions  within  the  CCS  (Chavez,  personal  communication). 

More  recently,  Belkin  et  al.  (2009)  carried  out  a  frontal  mapping  characterization  of  the 
CCS  based  on  available  SST.  Their  frontal  mapping  for  July  1985  and  December  1985  (see  their 
Eigure  14)  showed  a  dramatic  change  of  scale  (an  order  of  magnitude)  from  summer  to  winter. 
The  near-chaotic  (disorganized  and  almost  homogeneous)  distribution  of  fronts  in  the  summer 
changes  to  an  organized  long-filament  (-1000  km)  frontal  distribution  with  gaps  in  between  in 
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winter.  This  pattern  is  also  observed  in  other  eastern  boundary  current  regions  and  is  indicative 
of  the  prevalence  of  transitional  and  short-lived  fronts  associated  with  swift  dynamical  processes 
with  time-scales  of  1-10  days.  This  idea  supports  the  Mooers  and  Robinson  (1984)  construct  of 
meandering  jets  and  filaments  interacting  with  eddies  and  upwelling  regions  forced  by  short- 
period  atmospheric  winds.  A  more  recent  overview  of  the  processes  in  the  California  Current 
System  is  presented  by  Barth  and  Chekley  (2009). 

Modeling  and  assimilation  in  the  larger  California  Current  region  for  understanding  the 
complex  processes  have  been  a  major  focus  of  many  investigators,  including  McWilliams  et  al. 
(1985),  Miller  et  al.  (1999),  Mercheseilo  et  al.  (2003),  Penven  et  al.  (2006),  Caput  et  al.  (2008 
a,b),  Moore  et  al.  (2009)  and  others  listed  in  these  papers.  During  2003  and  2006,  major 
experiments  occurred  in  the  coastal  region  of  Monterey  Bay  (around  the  middle  of  our  coastal 
zone),  revealing  some  of  the  faster  upwelling  dynamics  processes  and  their  interactions  with 
larger-scale  features  including  offshore  eddies,  the  CUC  and  coastal  trapped  waves  (e.g.  Ramp  et 
al,  2008;  Chao  et  al,  2008;  Haley  et  al,  2009;  Shulman  et  al,  2009;  Ramp  et  al,  2010-  this-issue). 

It  is  clear  from  the  above  overview  that  this  region  is  dynamically  complex  and  rich  with 
short-term  energetic  phenomena.  Thus  it  is  important  to  (i)  develop  a  synoptic  description  of  the 
system  elements  and  (ii)  characterize  and  document  the  scales  and  variability  for  possible  use  in 
modeling,  initializing,  validating  and  assimilating  in  a  numerical  prediction  system. 

3.  Features,  Scales  and  Variability  in  the  CCS 

Based  on  the  above  overview,  we  identified  the  following  key  points  to  be  considered  for 

a  synoptic  synthesis  of  features:  (i)  a  synoptic  state  that  can  sustain  a  frontal  jet  of  higher 

velocity  similar  to  Brink  et  al.  (1991);  (ii)  the  low-salinity  core  of  the  CC  might  be  influenced  by 
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these  southward-flowing  synoptic  coastal  jets  mixing  with  the  mean  equatorward  flow  system 
(which  is  part  of  the  large-scale  subtropical  gyre);  (iii)  offshore  of  Monterey  Bay  at  around 
37. 5N,  the  synoptic  jet  (of  Strub  et  al.,  1991)  and  the  mean  current  (of  Collins  et  ah,  2003)  are 
within  one  degree  of  each  other;  (iv)  the  addition  of  a  number  of  eddies  might  provide  a  more 
stable  configuration. 

We  present  below  various  aspects  of  the  synoptic  variability  of  the  prevalent  circulation 
features  in  the  CCS,  focusing  on  the  summer  (July-August)  season.  Some  of  these  features  can 
be  seen  in  Figure  3  of  Strub  et  al.  (1991)  which  depicts  surface  pigment  concentration  from  the 
CZCS  satellite  data  from  June  15,  1981.  The  characteristic  CCS  surface  features  can  be  seen  in 
similar  SST  and  color  images  (Strub  and  James,  2000).  The  features  that  we  study  and 
characterize  are:  (i)  the  California  Current  (mean  and  its  core  jets);  (ii)  the  California 
Undercurrent;  (iii)  the  Inshore  Countercurrent  (Davidson  Current);  (iv)  the  coastal  transition 
zone  that  separates  the  shelf  circulation  from  the  offshore  flows;  (v)  the  eddies;  (vi)  the 
anomalous  pools;  (vii)  the  baroclinic  jets  along  upwelling  fronts;  (viii)  the  filaments;  (ix)  the 
squirts  and  (x)  the  mushroom  heads. 

In  particular,  we  provide  a  comprehensive  synthesis  of  the  typical  synoptic  widths, 
locations  (distance  from  the  coast),  vertical  extents,  and  core  characteristics  of  these  features. 
The  dominant  spatial  and  temporal  scales  of  variability  of  each  of  these  features  are  identified 
from  past  observational,  theoretical  and  modeling  analyses.  They  are  indicated  in  Table  1.  We 
discuss  some  of  these  features. 

3.1  The  California  Current 

The  current  understanding  is  that  the  core  of  the  CC  is  at  the  inshore  edge  of  the  current. 
This  low-salinity  core  is  associated  with  the  separating  coastal  jet  (Barth  et  al.,  2000)  at  the 
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surface.  This  frontal  jet  also  meanders  and  develops  eddies  as  it  progresses  southward  and 
interacts  with  the  onshore  upwelling  regions  and  their  filaments.  The  low  salinity  core  was 
probably  first  described  by  Lynn  and  Simpson  (1987).  They  defined  three  cross-shore  zones 
based  on  seasonal  signal  in  dynamic  height  field  from  the  CALCOFI  data.  Minimum  in  the 
dynamic  height  field  was  obtained  in  the  winter  in  the  offshore  region  due  to  cooling.  Maximum 
in  the  dynamic  height  was  observed  in  the  winter  caused  by  currents.  Eddies  dominated  the 
transition  zone  in  between,  where  the  standard  deviation  of  the  dynamic  height  was  high, 
although  its  seasonal  signal  was  relatively  small.  The  core  of  the  CC  lies  in  this  band.  This  is  also 
consistent  with  the  Strub  et  al.  (1991)  description  of  the  jet  flowing  through  the  field  of  eddies 
(See  Fig.  15  ofLS87). 

A  synergistic  set  of  three  simple  model  alternatives  for  the  flow  structure  of  cold 
filaments  proposed  by  Strub  et  al.  (1991)  (their  Figure  1)  is  shown  here  in  Figure  2.  They 
proposed  that  the  idealized  squirt  is  generated  by  nearshore  convergences,  such  as  those  caused 
by  local  wind  relaxations  around  capes,  so  the  squirts  have  an  onshore  energy  source. 
Alternatively,  the  field  of  mesoscale  eddies  could  draw  recently  upwelled  water  away  from  the 
coast,  which  is  an  offshore  energy  source.  The  third  alternative  was  to  imagine  a  continuous 
southward  jet  meandering  offshore  and  onshore.  This  jet  may  entrain  coastally  upwelled  water 
and  pull  it  offshore.  While  all  of  the  preceding  hypotheses  may  contribute  to  squirt  formation, 
the  last  hypothesis  of  the  meander-eddy-upwelling  is  becoming  the  prevailing  view  of  the  CCS 
in  spring/summer  (see  also  Ramp  et  ai,  2009,  Pringle  and  Dever,  2009). 

In  a  larger  context,  the  meander-eddy-upwelling  system  as  initially  proposed  by  Mooers 
and  Robinson  (1984)  and  adapted  here  in  Figure  3  as  a  synoptic  basis  for  synthesis  of  multiscale 
elements  should  be  able  to  reproduce  realistic  jets,  eddies  and  squirts  in  a  dynamical  modeling 
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framework.  The  schematic  in  Figure  3  shows  the  complexity  in  the  meander-eddy-upwelling 
system  of  the  California  Current  off  Oregon  and  California  in  the  spring/summer.  The  coastal  jet 
coming  from  the  north  moves  offshore  as  the  season  progresses  (denoted  by  the  March,  May, 
and  July  lines),  and  contributes  downstream  to  the  low-salinity  core  of  the  California  Current. 
This  core  delineates  the  higher  variability  region  nearshore  from  the  less  active  regions  offshore. 
The  offshore  region  consists  of  the  mean  southeastward  flow  of  the  wind-driven  subtropical 
gyre.  The  inshore  region  is  populated  by  upwelling  centers  concentrated  near  the  capes.  Eddies 
generated  by  the  meandering  jet  are  found  in  both  regions,  and  may  exhibit  either  cyclonic  or 
anticyclonic  rotation.  Temporally  transient  eddies  of  various  sizes  have  been  observed. 

3.2  The  California  Undercurrent  (CUC) 

The  availability  of  synoptic  data  for  the  CUC  is  sparse.  Differences  in  periods  of  analysis, 
geographical  coverage,  and  lack  of  high-resolution  synoptic  surveys  are  the  challenges  of 
pinpointing  the  location  and  behavior  of  this  feature.  Notably,  Collins  et  al.  (1996)  showed  that 
the  path  of  the  CUC  during  summer  1993  was  close  to  the  continental  slope  between  San 
Francisco  (37. 8N)  and  St.  George  Reef  (41. 8N).  Pierce  et  al.  (2000)  analyzed  in  detail  a  follow¬ 
up  survey  of  105  shipboard  ADCP  velocity  sections  during  July  to  August  1995.  They  found 
that,  on  monthly  average,  the  CUC  has  a  core  speed  of  >10  cm/sec  between  200-275  m  depth  at 
20-25  km  off  the  shelf  break.  These  sections  and  the  core  path  are  available  from  the  website 
http://diana.oce.orst.edu. 
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3.3  The  Coastal  Uywelling  Regions 


Seasonal  variations  in  the  wind  direction  provides  for  different  regimes  along  the  coast  of 
California.  From  March  through  September,  prevailing  northwesterly  winds  combined  with  the 
effect  of  the  earth’s  rotation,  drive  surface  waters  offshore,  leading  to  consistent  upwelling 
region  around  the  many  capes  (e.g.  Fratantoni  and  Haddock,  2009).  These  deep,  cold  nutrient- 
rich  waters  are  responsible  for  high  productivity  of  California’s  near-shore  region  (DSR  v53,  25- 
26,  2006).  The  upwelling  activity  subsides  in  September,  when  the  northwesterlies  die  down, 
which  leads  to  downwelling  in  October.  The  relaxation  of  the  upwelling  in  August  is  described 
by  Ramp  et  al.  (2005)  and  in  Rosenfeld  et  al.  (1994).  Haley  et  al.  (2009)  reported  four¬ 
dimensional  simulation  studies  of  dynamical  processes  in  the  Monterey  Region,  while  Ramp  et 
al.  (2010)  compare  different  simulations  of  such  dynamics,  pointing  to  the  importance  of  model 
resolution  and  of  longer  coastal  waves.  The  wintertime  change  in  the  atmospheric  pressure  leads 
to  southwesterlies  along  the  coast.  The  response  along  the  coast  is  generally  seen  as  a  northward 
flow  reversal  at  the  surface,  called  the  Inshore  countercurrent,  or  the  Davidson  Current. 

The  inshore  region  of  the  CTZ  is  complicated  by  the  poleward  flow  being  interrupted  by 
the  summer  upwelling  jets  flowing  equatorward  near  the  mouth  of  Monterey  Bay  (Collins  et  al., 
2003).  The  water-mass  characteristics  of  the  Monterey  Bay  during  the  upwelling  and  relaxation 
periods  were  described  in  detail  by  Warn-Vamas  et  al.  (2007).  For  the  development  of  the 
FORMS  for  the  CCS  region  during  summer,  we  focus  our  study  to  the  six  upwelling  regions 
around  the  six  capes:  Cape  Blanco,  Cape  Mendocino,  Pt.  Arena,  Pt.  Reyes,  Pt.  Ano  Nuevo  and 
Pt.  Sur. 
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3.4  Other  Mesoscale  Features  in  the  CCS  Region 


According  to  Collins  et  al.  (2003),  during  fall  and  winter,  when  the  upwelling  is  not 
present,  the  Inshore  Countercurrent  (Davidson  Current)  has  two  distinct  cores,  one  on  the  eoast 
and  one  50  km  offshore.  Sinee  our  focus  is  on  the  development  of  CCS-FORMS  for  the 
summertime,  we  defer  to  the  development  of  this  feature.  Apart  from  these  three  eurrents,  there 
are  very  energetie  and  important  shallow  water  features  in  the  CTZ  around  the  eoastal  region  of 
Monterey  Bay  and  the  California  Current  System.  The  filaments,  eoastal  eddies,  anomalous 
pools,  jets  along  the  upwelling  fronts  and  mushroom-like  vortiees  are  eritical  to  reproduee  for 
aeeurate  forecasting.  The  seales  and  variability  of  these  features  have  been  identified  in  Table  1, 
with  relevant  references  cited  in  Table  2. 

These  features  sueh  as  filaments,  jets,  and  mushroom-vortiees  are  related  to  specifie 
dynamical  processes,  ineluding  instability,  upwelling  and  wind-driving.  It  has  been  proposed  by 
Brink  et  al.  (1991)  that  subduetion  plays  a  very  important  role  inside  the  filaments,  foreing  these 
cold  filaments  downwards,  below  the  ambient  warm  surface  waters.  The  mushroom  vortices  are 
generated  only  if  certain  eonditions  of  the  flow  field  are  met  (Mied  et  al.,  1991). 

For  this  first  implementation  of  FORMS  for  the  CCS,  we  have  exeluded  sueh 
eomplieations.  However,  as  will  be  shown  later,  onee  the  primary  features  are  in  plaee,  the 
ensuing  model  integration  generates  mesoscale  and  submesoscale  features  such  as  filaments  and 
mushroom  vortices.  Examples  of  such  realistic  simulations  were  given  recently  for  the 
southwestern  Atlantie  region  by  Calado  et  al.  (2008). 
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4.  Feature  Models  for  Different  Features  in  the  CCS 


A  prerequisite  synthesis  of  the  prevalent  features  in  the  CCS  was  presented  in  Section  3 
and  their  spatial  and  temporal  scales  are  summarized  in  Table  1.  In  particular,  we  have 
investigated  the  typical  synoptic  width,  location  (distance  from  the  coast),  vertical  extent,  and 
core  characteristics  of  these  features.  For  the  CCS-FORMS,  we  restrict  ourselves  to  a 
combination  of  the  CC,  CUC  and  upwelling  feature  models  in  a  climatological  background  for 
initialization.  Eddies  can  be  added  in  a  synoptic  forecasting  system,  and  multiple  examples  of 
such  eddy  feature  models  can  be  found  in  the  studies  by  Gangopadhyay  and  Robinson  (2002), 
Calado  et  al.  (2006),  Calado  et  al.  (2008)  and  others.  In  this  section  we  present  the  feature 
models  that  were  developed  based  on  previous  synoptic  studies  for  a  subset  of  those  features  for 
this  first  application  to  numerical  modeling. 

4.1  Feature  Models  for  the  California  Current  and  its  Low  Salinity  Core 

The  schematic  in  Figure  3  shows  the  setup  of  different  features  in  the  CCS  region.  From 
a  synoptic  viewpoint,  the  California  Current  mean  flow  has  a  meandering  jet  structure  at  its  core. 
We  chose  the  low-salinity  pool  in  the  core  of  the  CC  as  the  primary  characteristic  and  developed 
a  feature  model  for  this  structure.  Low  salinity  is  a  reasonable  descriptor  of  the  path  of  the  CC 
(see  Figures  9  and  1 1  of  Lynn  and  Simpson,  1987).  The  salinity  minimum  is  offshore  of  the 
velocity  maximum.  Lynn  and  Simpson  (1987)  suggested  that  this  is  due  to  more  mixing  with 
inshore  waters  than  with  offshore  waters. 

Ligure  4  shows  the  CalCOLI  station  locations  and  the  vertical  sections  of  mean  salinity 
for  July  along  CalCOLI  stations  60,  90  and  1 10.  The  low-salinity  pool  is  visible  in  these  sections 
as  a  distinct  core.  We  choose  to  model  the  isohalines  in  the  low-salinity  core  using  a  series  of 
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parabolas  having  their  foci  on  the  same  vertical  (see  Figure  5a).  This  vertical  line  passes  through 
the  center  of  the  core,  is  the  axis  of  all  the  parabolic  isohalines  and  is  perpendicular  to  the 
directrix  of  the  set  of  parabolas. 

The  parabolic  formulation  for  the  low  salinity  pool  for  a  point  (x,  z)  (Fig.  5a)  is  given  by: 
x^=Aa{z  +  H)  (1) 

where  H  is  the  depth  of  the  isohaline  core,  and  the  parabolic  isohaline  has  vertex  at  (0,  -H)  in  the 
X-Z  plane  (Fig.  5a). 

At  the  surface,  where  z  =  0,  we  can  use  the  salinity  data  obtained  by  Lynn  et  al.  (1982) 
and  presented  by  (Lynn  and  Simpson,  1987)  for  evaluating  the  “a”  values  of  the  parabolic 
isohalines. 

Specifically,  for  the  surface  points  of  each  parabolic  isohaline,  one  can  write  a  simple 
relationship  between  the  maximum  width  (W)  of  the  pool  and  the  depth  {H)  of  the  isohaline  as: 
(W72)^  =  4a//;  which  leads  to  a  =  (!///)( W/4)^.  Using  the  observed  width  and  depth  of  each  of  the 
isohalines  for  each  of  the  sections  in  Figure  4,  the  “a”  values  of  these  parabolas  were  obtained. 
These  data  have  been  tabulated  in  Table  3.  The  modeled  parabolic  isohalines  corresponding  to 
Line  60  are  shown  in  Figure  5b,  and  they  compare  reasonably  well  with  the  observed  isohalines 
shown  in  Figure  4b. 

We  next  distinguish  between  the  “core”  area  of  the  pool  and  the  parabolic  isohalines 
outside  the  core.  The  salinity  in  the  core  area  is  assumed  to  be  constant.  So,  from  surface  to  the 
core  depth  (Co)  bounded  by  the  parabolic  isohaline  of  So  in  the  x-z  vertical  section.  Thus,  from  a 
feature  model  perspective,  given  a  specific  set  of  So,  Co,  we  now  need  to  find  the  salinity  at  each 
point  outside  the  “core”  within  the  pool.  Equation  (1)  above  defines  boundary  of  the  core. 
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As  we  have  seen  above,  the  parabolic  assumption  allows  one  to  compute  the  “a”  value  for 
different  “/f  ’  and  once  these  are  plotted,  it  was  noted  that,  for  increasing  value  of  H,  “a” 
increases  almost  linearly,  which  led  us  to  relate  and  ‘a’  by  the  following  equation 


H  = 


\  da  j 


(a)  +  Cq 


(2) 


where  Co  is  the  core  depth. 

Plugging  the  value  of  H  in  equation  1  we  have 


xf  =  4a 


Zi  + 


^  dH^ 


\  da  j 


(a)  +  Co 


Rearranging  terms,  a  quadratic  in  “a”  can  be  formed 

^  dH^ 


\  da  j 


(a  )  +  (Zi  +Co)(a)- 


f  2  A 
v4y 


=  0 


whose  real  solution  is  given  by: 


a  = 


r  ,  f dH^  ,  ~ 

1/2 

-  (Zj  -1-  Cq)  -1- 

(z,+Q)^+  —  (vf) 

\  uci  J 

^  dH^ 
\  da  J 


(3) 


Thus,  once  we  get  an  ‘a’,  we  can  use  equation  2  to  find  out  the  depth,  H,  of  the  isohaline.  In  this 
way  the  isohaline  gets  defined.  Analysis  of  the  observed  core  (from  Figure  4)  at  several  lines 

(  dS"\ 

yielded  an  almost  constant  value  of  —  .  Thus,  similar  to  H  of  an  isohaline,  point-wise  salinity 

\da  ) 

can  now  be  obtained  using  the  equation 


S  = 


\da  J 


(a)  +  Sq 


(4) 
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where  So  is  the  salinity  of  the  core  of  the  pool.  So,  it  turns  out  that  we  can  express  S  as  a 

function  of  So,  Co,  —  and  -  .  Given  these  four  parameters,  we  can  now  determine  the 

da  da 

salinity  (using  Table  3,  Equations  3  and  4)  at  any  point  in  the  low-salinity  pool. 

Figure  6a-b  shows  the  temperature  and  salinity  profiles  for  this  core,  which  were  adapted 
from  Lynn  and  Simpson’s  (1987)  analysis  of  the  data  obtained  through  the  long-term  monitoring 
program  CalCOFI .  From  such  a  typical  salinity  profile,  one  can  extract  the  parameters  So  and  Co 
for  the  feature  model.  The  salinity  and  depth  variation  with  “a”  are  given  by  the  numbers  in 
Table  3.  For  the  implemented  feature  model,  the  initial  path  of  the  CC  is  a  representative  one, 
and  would  have  to  be  obtained  on  the  basis  of  SST  or  color  signatures.  It  is  suggested  that  the 
core  of  the  CC  is  better  identified  from  color  images  (such  as  those  from  MODIS)  as  the 
chlorophyll  boundary.  In  the  future,  when  Sea  surface  salinity  fields  will  be  available  from 
NASA’s  Aquarius  mission,  one  would  expect  to  identify  this  core  of  the  CCS. 

4.2  Feature  Models  for  the  California  Undercurrent  ( CUC) 

The  scarcity  of  data  for  the  CUC  was  highlighted  in  the  overview  section.  The  only 
complete  and  continuous  profiling  was  done  by  Pierce  et  al.  (2000),  which  is  shown  in  Figure  7. 
It  is,  however,  worthwhile  to  note  that  the  location  of  the  CUC  offshore  of  the  Monterey  Bay 
region  (the  focus  of  AOSN-II  during  summer  2003)  was  highly  variable.  Early  observations  of 
the  CUC  in  this  region  (36N  to  37N)  by  Wickham  (1975)  indicated  a  complex  alongshore  flow 
near  the  coast,  in  both  poleward  and  equatorward  directions  in  the  upper  500  m  offshore  close  to 
123W.  In  particular,  he  observed  a  poleward  flow  being  subdivided  into  two  parts  by  a  strong 
(60  cm/sec)  equatorward  jet  down  to  500  meters  (Fig.  11  of  Wickham,  1975).  The  recent  NMFS 
survey  sections  (Fig.  3)  by  Pierce  et  al.  (2000)  indicate  a  maximum  of  the  CUC  shifted  offshore 
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at  36.47N  and  at  36. 8N,  while  flowing  very  close  to  the  coast  north  of  37N  (see  Fig.  4b  of  Pierce 
et  al.,  2000).  A  section  at  35.97N  shows  a  clear  equatorward  jet  in  the  middle  of  the  weak 
poleward  flow. 

In  fact,  Rischmiller  (1993)  found  that  the  area  of  poleward  flow  off  of  Point  Sur  extended 
well  beyond  the  upper  continental  slope  up  to  a  distance  of  200  km  from  the  coast  and  to  a  depth 
of  750-1000  m  or  greater  (Garfield  et  al.,  1999).  Based  on  Lagrangian  drifters  during  1992-95, 
Garfield  et  al.  (1999)  plotted  the  ensemble  and  individual  spaghetti  diagrams  (Figs.  5  and  6  of 
their  paper),  which  clearly  show  that  the  CUC  flow  is  further  offshore  across  the  Monterey  Bay 
region  than  it  is  to  the  north  of  this  region. 

With  this  background,  we  have  identified  the  July  1984  CalCOFI  data  set  (analyzed  by 
Pierce  et  al.  2000)  as  a  candidate  analysis  data  set  for  the  basis  of  a  CUC  template.  A  first-order 
feature  model  for  the  CUC  will  be  implemented  on  the  basis  of  this  data  set  from  July  1984,  and 
other  supporting  temperature-salinity  observations.  These  feature  model  profiles  are  shown  in 
Figure  8.  Preliminary  analysis  suggests  that  a  warm  and  saline  water  mass  can  be  traced  at  200- 
300  m  depth.  Curiously,  the  path  of  this  water  mass  follows  the  2500  m  isobath  in  the  latitudinal 
range  of  34-37. 5N.  This  offshore  trajectory  resembles  the  single  drifter  trajectory  of  Garfield  et 
al.  (1999),  as  was  shown  in  their  Figure  6. 

At  depth,  the  CUC  is  observed  to  shed  eddies.  Garfield  et  al.  (1999,  2001)  tracked  38 
approximately  isobaric  RAFOS  floats  ballasted  at  300m  in  the  CUC.  They  found  2  size-classes 
of  eddies:  small  (radius  <  35  km)  rapidly  anticyclonically  rotating  submesoscale  coherent 
vortices  (that  they  called  “cuddies”),  and  larger  (radius  >50  km)  eddies  formed  due  to  other 
processes  such  as  baroclinic  instability.  They  found  that  the  cuddies  propagate  westward  into  the 


19 


interior  with  speeds  of  1-2  cm/s  (see  their  Figure  5).  We  will  show  later  that  the  feature-modeled 
cue  develop  eddies  similar  to  these  at  subsurface  depths  during  the  dynamical  simulations. 

4.3  Coastal  Upwelling  Front  Feature  Model 

This  subsection  documents  the  feature  models  of  the  coastal  upwelling  fronts.  Six 
locations  are  chosen  for  the  upwelling  sites  -  Cape  Blanco,  Cape  Mendocino,  Pt.  Arena,  Pt. 
Reyes,  Pt.  Ano  Nuevo  and  Pt.  Sur.  For  the  estimate  of  the  upwelling  front  feature  model,  two 
temperature  profiles  are  required:  one  for  inshore  and  the  other  for  offshore.  Base  data  were 
obtained  from  various  sources  and  various  cruises  listed  in  Table  1  of  Kim  et  al.  (2007). 

The  upwelling  feature  models  for  each  cape  follow  the  methodology  of  Shaji  and 
Gangopadhyay  (2007).  An  upwelling  feature  model  requires  the  specification  of  two  T  and  S 
profiles  -  an  inshore  profile  and  an  offshore  profile,  which  are  dynamically  interpolated  with  a 
tank  function  (see  Shaji  and  Gangopadhyay  (2007)  and  Calado  et  al.  (2008)  for  details). 

The  upwelling  region’s  temperature  is  given  by, 

T(q;  z)  =  To(z)  -r  [Tu(z)  -  To(z)]  m(q,  z)  (5) 

where 

m(q,  z)  =  Vi  +  Vi  tanh  [(q  -  0z)/A.] 

is  a  melding  function,  rj  is  the  cross-frontal  distance  from  the  axis  of  the  front,  and  z  is  positive 
vertically  upward.  T{rj,  z)  is  the  upwelling  regional  temperature  distribution,  Tu(z)  is  the 
upwelling  nearshore  temperature  profile,  and  To(z)  is  the  offshore  temperature  profile.  6^  is  the 
slope  of  the  front,  and  X  is  the  e-folding  half- width  of  the  front  (=  L/2)  (see  Figure  9). 

Figure  10  shows  the  coastal  upwelling  front  FM  T  profiles  (left-  Fig.  10a)  and  newly 
generated  S  profiles  (Fig.  10b)  at  Cape  Blanco.  Figure  10c  is  a  TIS  diagram  where  inshore  and 
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offshore  T  and  S  (red)  and  derived  FM  T/S  profiles  (blue)  are  compared  to  show  the  range 
captured  by  the  feature  model.  The  resulting  upwelling  section  at  initialization  is  shown  in 
Figure  11.  Clearly  such  initialization  would  have  significant  impact  on  the  predictive  capability 
of  a  forecast  system  in  the  shallow  coastal  region  of  the  CCS.  Similar  constructs  of  the  upwelling 
FMs  were  derived  and  implemented  for  the  other  five  Capes. 

5.  Synoptic  Initialization,  Model  Domain  and  Numerical  Parameters 

As  indicated  in  Section  1,  feature-oriented  methodology  requires  the  development  of  a 
synoptic  circulation  template  for  the  region.  This  template  is  designed  from  a  feature-based 
synthesis  of  the  regional  circulation  pattern.  From  it,  a  map  of  strategic  sampling  locations  for 
placing  feature  model  profiles  is  produced;  it  provides  "synthetic  synoptic  expressions"  for 
fronts,  eddies,  gyres,  and  other  circulation  structures  and  water  masses  at  the  initialization  or 
model-updating  phase.  A  suitable  climatology  is  then  selected  to  provide  the  background  of  the 
synoptic  features.  Further,  a  multiscale  objective  analysis  (OA;  Lermusiaux,  1999a,b)  is 
performed  to  meld  the  synoptic  and  large-scale  fields,  which  results  in  a  three-dimensional  field 
for  initializing  a  numerical  model  grid.  In  this  section,  we  describe  the  initialization  procedure, 
followed  by  the  details  of  the  selected  model  domain  and  its  parameters. 

5.1.  Construction  of  Synoptic  Initialization  Using  the  Feature-based  Multi-scale 
Synthesis 

Figure  12  shows  an  example  of  the  synoptic  circulation  template  for  the  CCS  region.  The 
major  advantage  in  using  such  a  template  is  the  ability  to  resolve  synoptic  structures  on  the  basis 
of  past  oceanographic  knowledge  of  the  region,  even  when  there  is  a  lack  of  observations  for  the 
synoptic  features.  The  red,  black  and  blue  stars  in  Figure  12  represent  the  high-resolution 
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sampled  stations  where  the  feature  model  (FM)  temperature  and  salinity  profiles  are  provided  to 
resolve  three  individual  synoptic  features.  In  this  region,  the  CCS  (red  stars),  the  CUC  (black 
stars)  and  the  upwelling  regions  (blue  stars)  are  the  primary  synoptic  features  that  are  represented 
in  the  circulation  template.  The  axis  of  the  core  of  the  CC  is  determined  as  a  line  parallel  to  a 
smoother  version  of  the  coastline  at  about  500km  offshore.  The  locations  of  the  perimeter  of  the 
upwelling  regions  are  all  based  on  historical  observations  around  the  capes.  Note  that  only  a 
limited  number  (-100)  of  stations  are  needed  to  describe  the  synoptic  behavior  of  this  region. 
These  are  again  placed  based  on  a  handful  of  FM  parameters  described  in  Section  4. 

Once  the  feature  models  are  developed  (as  described  in  Section  4),  an  initialization  field 
is  then  generated  by  objectively  analyzing  the  synthetic  FMs  on  the  synoptic  template  in  a 
selected  background  climatology  for  a  multi-scale  synthesis.  The  green  stars  in  the  circulation 
template  in  Figure  12  are  the  locations  of  the  Boyer  et  al.  (2005)  14  -degree  climatology  data. 

The  multiscale  OA  is,  in  fact,  done  in  two  stages.  In  the  first  stage,  the  largest  dynamical  scales 
are  resolved  at  each  level,  using  estimated  large-scale  e-folding  spatial  decays,  zero-crossings 
and  temporal  decay.  The  background  for  the  first  stage  OA  is  the  horizontal  average  of  the 
climatological  data.  In  the  second  stage  of  the  OA,  the  synoptic  mesoscale  phenomena  are 
resolved  using  their  estimated  space-time  scales.  The  background  for  this  OA  is  the  first  stage 
OA.  Once  the  initialization  field  is  available  as  a  result  of  the  multi-scale  synthesis,  any 
numerical  model  can  then  be  used  for  nowcasting,  forecasting  and  dynamical  process  studies 
with  such  initialization  fields. 

We  use  Boyer  et  al.'s  (2005)  temperature  and  salinity  fields  objectively  analyzed  for 
August  with  14  degree  of  spatial  resolution.  The  method  of  objectively  analyzing  the  data  is  the 
same  as  outlined  in  Boyer  et  al.  (2005)  for  temperature  and  salinity.  There  are  differences  in 
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technique  necessitated  by  the  different  scales  of  the  analysis.  In  the  one-degree  analysis,  the 
radius  of  influence  is  250  km;  in  the  quarter-degree  analysis,  the  radius  of  influence  is  40  km. 

5.2.  The  POM  Parameters  and  Domain  Setup 

Initialization  of  any  dynamical  model  requires  the  specification  of  each  prognostic 
variable  at  each  grid  point  of  the  three-dimensional  model  domain.  In  this  work  we  chose  the 
terrain-following  coordinate  system  for  this  implementation  of  the  Princeton  Ocean  Model  - 
POM  (Ezer  and  Mellor,  1994;  Blumberg  and  Mellor,  1987  and  Ezer  and  Mellor,  1997).  The 
domain  consists  of  150  x  150  grid  points  in  the  horizontal  with  8  km  resolution  and  31  vertical 
levels.  The  model  is  limited  geographically  by  34.38N,  48.82N  and  -138. 22E,  -118.68E. 

The  model  bathymetry  is  based  on  the  ET0P05  data  set.  The  depth  range  is  5  to  6000  m. 
On  the  north,  south  and  west  we  used  open  boundary  conditions.  Along  these  boundaries  the 
model's  thermohaline  structure  was  held  to  initial  temperature  and  salinity  values  via  partially- 
clamped  boundary  conditions  (Blumberg  and  Kantha,  1985).  Also,  a  radiation  condition  based 
on  Summerfield's  equation  was  applied  to  the  baroclinic  velocities  on  the  open  boundaries. 
Elather  radiation  condition  (Palma  and  Matano,  1998)  was  applied  to  the  barotropic  velocities, 
and  a  non-gradient  condition  for  the  free  sea  surface  elevation  was  employed  as  well. 

As  a  consequence  of  the  relaxation  in  the  temperature-salinity  fields,  we  imposed  a 
geostrophic  flow  at  the  model's  open  boundaries  through  the  inclination  of  the  isopycnal 
surfaces.  In  other  words,  the  climatological  signal,  rendered  by  the  inclinations  of  the  isothermal 
and  isohaline  surfaces,  contains  the  signal  of  both  the  velocity  component  generated  by  wind 
forcing  and  the  component  associated  with  the  thermohaline  circulation.  The  wind  forcing  was 
generated  from  the  mean  August  field  extracted  from  QuikScat  data  (Graf  et  al.,  1988).  We  also 
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included  the  deep  CUC,  underlying  the  CC,  in  the  regional  modeling,  as  well  as  the  influence  of 
the  thermohaline  circulation  component  on  the  mesoscale  activity  in  the  region  of  interest. 

5.3  Initialization  Strategy  for  Short-term  Forecasting  with  POM 

Synoptic  short-term  simulations  require  three-dimensional  specifications  of  initial  fields 
of  T,  S,  u  and  v.  However,  the  standard  practice  for  basin-scale  numerical  modeling  exercises  is 
to  specify  T  and  S  only,  with  zero  velocity  (Ezer  and  Mellor,  1994;  Ezer  and  Mellor,  1997; 
Mathew  and  McClean,  2005).  The  basin-scale  approach  requires  models  to  develop  the  inertia 
fields  in  a  one-to-ten  years  long  "spin-up"  period.  As  mentioned  earlier,  our  goal  is  to  develop  a 
system  for  short-term  synoptic  forecasting  where  we  cannot  afford  to  run  models  over  a  long 
spin-up  period.  Thus,  specifying  the  best  possible  velocity  field  in  addition  to  the  mass  (T  and  S) 
field  becomes  mandatory  for  synoptic  short-term  simulations.  Adjustments  to  wind  forcing  and 
baroclinic  internal  adjustment  are  also  necessary  for  synoptic  forecasting. 

Previous  work  on  such  adjustment  procedures  included  the  Harvard  Ocean  Prediction 
System  (HOPS)  in  the  studies  of  Gangopadhyay  et  al.  (1997,  2003)  and  Brown  et  al.  (2007  a,  b), 
the  Regional  Ocean  Modeling  System  (ROMS)  by  Shaji  and  Gangopadhyay  (2007)  and  Kim  et 
al.  (2007),  and  a  first  application  with  POM  by  Calado  et  al.  (2008).  The  strategy  for  setting  up 
POM  in  the  CCS  region  was  carried  out  using  climatological  temperature  and  salinity  initial 
fields  following  Ezer  and  Mellor’ s  (1997)  methodology,  with  the  additional  implementation  of 
geostrophic  velocity  as  initial  condition  for  the  baroclinic  velocities,  balanced  with  initial  T-S 
field.  Specifically,  we  held  the  T  and  S  fields  fixed  over  the  first  five  days  of  simulation  in  the 
diagnostic  mode,  an  approach  described  by  Ezer  and  Mellor  (1997).  These  authors  used  this 
diagnostic  mode,  available  in  POM,  to  dynamically  adjust  the  momentum  field  with  the 
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underlying  mass  field  and  to  the  topography.  A  similar  technique  for  initializing  POM  with 
geostrophic  velocities  was  used  by  Castelao  et  al.  (2004)  to  study  shelf  break  upwelling  in  the 
region  of  southeastern  Brazil.  The  initial  geostrophic  velocity  was  computed  from  the 
temperature  and  salinity  by  the  OA  methodology  described  by  Lozano  et  al.  (1996).  A 
representative  initial  field  of  temperature,  salinity  and  velocity  can  be  seen  for  day  1  in  Figure  13 
(left-top)  and  Figure  14  (left- top). 

6.  Dynamical  Simulations  with  Feature  Models 

The  long-term  goal  of  this  study  is  to  develop  the  capability  of  nowcasting  and 
forecasting  using  the  FORMS  technique  in  the  CCS  region.  So  far,  we  have  discussed  the 
development  of  a  set  of  appropriate  feature  models  for  this  region  and  their  incorporation  into  a 
multiscale  objectively  analyzed  initialization.  In  this  section  we  describe  the  results  of  a  40-day 
simulation  with  such  initialization  using  the  POM  setup  described  above. 

6.1  The  Surface  Temperature  and  Salinity  Evolution 

Figures  13  and  14  show  the  evolution  of  the  CCS-FORMS  initialization  for  a  typical 
summer  synoptic  realization.  Temperature  fields  (with  current  vectors  superposed)  for  days  1,  5, 
11,  21,  31  and  41  are  shown  in  the  six  panels  in  Fig.  13  as  marked.  The  initial  field  shows  the 
upwelling  around  the  capes  and  the  southward-flowing  current  with  the  low-salinity  pool  (Figure 
14a)  identified  as  the  CC  core.  A  continuous  southward-flowing  current  develops  along  the  coast 
by  day  1 1 .  Additionally,  a  narrow  region  similar  to  the  coastal  transition  zone  develops  during 
the  third  week. 

As  mentioned  in  the  previous  section,  the  upwelling  regions  were  included  at 

initialization  as  feature  models.  Interestingly,  there  are  two  mesoscale  eddy  development  events 
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at  42N  and  at  38N.  Based  on  CTZ  satellite  and  field  data  analyses  during  May-July  1987  and 
1988,  Strub  et  al.  (1991)  determined  the  existence  of  a  surface  meandering  jet  surrounded  by  a 
field  of  eddies  between  36  and  42N.  Large  dynamic  height  ranges  corresponding  to  more  well- 
developed  eddies  were  found.  The  authors  maintained  that  the  eddies  occurred  mostly  inshore  of 
the  jet  north  of  39N,  and  were  both  inshore  and  offshore  of  the  jet  south  of  39N.  In  our 
simulations,  starting  with  the  jet  and  upwelling  regions  in  the  FORMS  setting,  a  large  eddy 
inshore  of  the  jet  is  formed  around  42N.  A  number  of  additional  eddies  (both  cyclonic  and  anti- 
cyclonic)  were  also  generated  during  the  simulation  (Fig.  13d,  e,  f).  On  the  larger  scale,  it  is 
interesting  to  note  that  the  dynamical  evolution  of  the  temperature  field  is  similar  to  the  synoptic 
SST  shown  in  Figure  la,  indicating  the  robustness  of  the  model  simulations. 

The  salinity  evolution  is  shown  in  Figure  14.  Note  how  the  low-salinity  core  of  the  CCS 
is  initialized  as  a  jet  feature  model  and  how  it  dynamically  adjusts  and  evolves  during  the 
simulation  while  interacting  with  the  surrounding  waters.  The  high  salinity  signature  of  the  eddy 
development  at  42N  is  clearly  seen  within  the  spread  of  the  low-salinity  coastal  waters.  The 
boundary  of  the  coastal  transition  zone  is  probably  best  designated  by  the  dynamic  narrow  band 
of  high-salinity  between  the  low-salinity  core  of  the  CC  and  upwelling  regions  along  the  coast. 
The  mesoscale  eddy  evolution  during  the  first  month  of  simulation  at  the  surface  is  presented  in 
Figure  15  (a-f)  in  terms  of  eddy  kinetic  energy.  Figure  16  shows  the  evolution  of  this  mesoscale 
field  in  a  near  coastal  focused  region. 

Several  realistic  features  are  observed  during  the  5-30  days  of  simulation  as  presented  in 
Figure  16.  First,  an  anticyclone  off  Monterey  Bay  is  seen  by  day  5  at  37N.  This  is  accompanied 
by  a  small  paired  cyclone  at  38N.  The  cyclone  continues  to  grow  and  maintain  itself  through  day 
30.  The  anticyclone  off  Monterey  Bay  interacts  with  the  meandering  jet  and  develops  a  three- 
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part  anticyclone-cyclone-anticyclone  system  by  day  20,  which  persists  through  day  30.  In  the 
northern  part  of  the  domain,  a  cyclone  at  42-43N,  1 26-1 27 W  is  apparent  by  day  10,  which  moves 
NW  continuously  through  day  30.  This  eddy  is  originally  connected  with  the  upwelling  at  Cape 
Mendocino  (41. 5N).  A  quick  succession  of  cyclones  and  anticyclones  are  observed  near  40N. 
This  region  near  40N,  126W  can  be  identified  as  the  most  variable  region  in  the  simulation. 

Specifically,  the  series  of  cyclones  and  anticyclones  developed  by  day  20  (Figure  16,  left- 
bottom)  warrants  further  discussion.  This  evolution  resembles  the  schematic  in  Figure  2:  the 
middle  as  well  as  the  3'^'^  panel  indicates  that  the  dynamical  evolution  of  the  CCS  is  a 
combination  of  the  series  of  eddies  and  the  jet-meander-eddy  system.  The  simulated  convergent 
shelf  flow  near  Cape  Mendocino  was  observed  and  described  by  Magnell  et  al.  (2003).  An 
anticyclonic  eddy  similar  (albeit  smaller)  to  the  observed  seasonally  recurring  giant  eddy 
recently  described  by  Crawford  et  al  (2009)  is  simulated  near  39N,  125W. 

Next  section  describes  the  evolution  of  the  undercurrent  and  its  associated  eddies  which 
are  shown  in  Figures  17  through  19.  Figure  20  shows  a  Chlorophyll  image  (for  September  2004) 
along  with  the  model  evolution  on  day  20.  Interestingly,  there  are  distinct  resemblances  in  the 
eddy  activity  in  a  number  of  sub-regions  between  the  independent  image  and  the  simulated 
temperature  field  which  was  evolved  by  POM  starting  from  the  feature-based  initialization 
without  any  eddy  field  (See  Figure  caption  for  details).  Such  comparisons  will  be  quantified  in  a 
hindcast  mode  in  future  for  multiple  hindcasts  to  help  develop  an  operational  forecasting  system 
for  the  CCS. 

6.2  Subsurface  Evolution  (300  m  and  1000  m) 

At  300  m,  the  known  circulation  of  the  CCS  during  summer  is  primarily  eddy-driven, 

with  significant  subsurface  eddies  (Brink  et  al.,  2000).  Figure  17  shows  the  evolution  of  the 
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temperature  salinity  and  velocity  fields  at  300  m.  The  subsurface  eddies  are  visible  in  both 
temperature  and  salinity  fields.  Since  the  low-salinity  pool  of  the  CC  core  was  feature-modeled 
to  be  restricted  vertically  above  160  m,  it  does  not  penetrate  to  the  intermediate  depths.  The 
subsurface  signature  of  the  surface  eddies  at  42N  and  38N  is  visible  in  Figs.  17a-b.  Figure  18 
shows  the  mesoscale  evolution  near  the  coast  inshore  of  DOW  during  the  first  four  weeks  of  the 
simulation.  The  subsurface  signature  of  the  surface  eddies  at  36N  (the  ACE-CE-ACE  triad)  is 
evident  on  day  20.  By  day  30  this  evolves  into  a  CE-ACE  pair  similar  to  that  at  the  surface. 

The  1000  m  flow  evolution  is  presented  in  Eigure  19  for  days  1,  5,  10,  20,  30  and  40.  It  is 
clear  from  the  simulation  that,  in  addition  to  forming  subsurface  eddies,  there  is  a  natural 
tendency  of  poleward  flow  development  along  the  California  coast.  Note  that  this  coastal 
undercurrent  is  the  natural  pathway  for  the  low-oxygen,  nutrient-rich  water,  which  originates  in 
the  eastern  tropical  Pacific.  Thus,  the  EORMS-CCS  could  be  useful  for  predicting  the 
biogeochemistry  of  the  California  Current  System  in  the  future. 

A  series  of  anticyclonic  eddies  developed  near  the  coast  inshore  of  the  poleward  flow 
during  the  simulation.  By  day  5,  an  ACE  develops  first  near  the  coast  at  37N.  A  poleward 
continuous  flow  develops  at  1000m  by  day  10.  The  ACE  at  37N  grows  between  day  10  and  day 
40.  By  day  20,  a  second  ACE  is  formed  inshore  of  poleward  flow  at  around  40-41N,  which 
grows  in  size  by  day  30.  A  third  inshore  ACE  is  seen  to  develop  at  44N  around  day  40.  In 
addition  to  the  ACEs,  an  ACE-CE  pair  was  observed  between  36  and  38N  (similar  to  the  surface 
pair)  on  day  20.  The  smaller  CE  signature  dissipated  by  day  30.  Table  4  summarizes  the  sizes 
and  nature  of  the  simulated  eddies  discussed  above  along  with  some  observational  evidence  of 
these  eddies  in  past  studies. 
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7.  Conclusion 


This  study  described  the  development  of  a  first-order  feature-oriented  regional  modeling 
system  (FORMS)  for  the  California  Current  system.  The  synoptic  circulation  features  include  the 
mean  flow  and  southeastward  meandering  jet(s)  of  the  California  Current,  the  poleward 
California  Undercurrent  and  six  upwelling  regions  along  the  western  coast  of  the  U.S.  (Figure  3). 
We  identified  the  typical  synoptic  width,  location,  vertical  extent,  and  core  characteristics  of 
these  features  and  their  dominant  spatial  and  temporal  scales  of  variability  from  past 
observational,  theoretical  and  modeling  studies,  which  are  listed  in  Table  1.  These  synoptic 
features  were  then  melded  with  available  climatology  of  the  region  for  initialization. 

Dynamical  simulations  over  a  40-day  period  using  POM  illustrate  the  applicability  of 
FORMS  to  a  transient  eastern  boundary  current  region  such  as  the  CCS.  The  simulations 
initialized  with  FORMS  realistically  evolved  various  interactions  among  the  low  salinity  pool, 
the  background  climatology  and  the  upwelling  regions.  A  number  of  anticyclones  and  cyclones 
appeared  at  places  where  observational  evidences  of  such  eddies  were  documented  in  previous 
studies.  These  are  listed  in  Table  4.  In  addition,  a  countercurrent  develops  in  the  intermediate  to 
deep  waters  along  the  coast,  which  sheds  realistic  eddies  (“cuddies”)  near  a  number  of  capes. 

This  physical  phenomenological  validation  is  the  first  step  towards  building  a  FORMS- 
based  prediction  system  for  this  region.  The  next  step  would  be  to  statistically  validate  such 
simulations  (initialized  with  low-salinity  pool,  coastal  upwelling,  CUC,  CC,  eddies  and 
background  climatology)  following  the  works  of  Robinson  and  Gangopadhyay  (1997)  for  the 
western  North  Atlantic.  Calibration  of  such  a  system  for  MB03  and  MB06  and  other  hindcast 
studies  will  lead  to  a  robust  implementation  of  an  operational  system  for  the  CCS  in  the  near 
future. 


29 


Acknowledgement 


This  work  was  funded  by  the  Office  of  Naval  Research  grants  NOOO 14-03 -1-04 11  and  N00014- 
03-1-0206  at  the  University  of  Massachusetts  at  Dartmouth.  Leslie  Rosenfeld’s  participation  was 
supported  by  ONR  grant  N00014-03-WR- 20009.  Our  special  thanks  to  Dr.  Jim  Bellingham  and 
his  group  at  MBARI  for  numerous  informative  discussions  and  making  the  data  available  from 
the  MUSE  and  AOSN2  experiments.  We  are  grateful  to  all  who  contributed  to  this  effort  by 
providing  datasets.  We  would,  especially  like  to  acknowledge  Drs.  Ken  Brink,  Tommy  Dickey, 
Francisco  Chavez,  Toby  Garfield,  Alex  Warn-Vamas,  Shaun  Johnston,  Curt  Collins,  Margaret 
McManus,  Libe  Washburn,  and  Erika  McPhee-Shaw  for  their  help  and  support  in  procuring  the 
datasets  from  various  sources.  We  appreciate  the  editorial  assistance  provided  by  Mr.  Frank 
Smith  and  graphics  help  of  Ms  Carolina  Nobre  at  SMAST.  This  is  our  final  AOSN  Contribution. 


30 


References 


Abbott,  M.R.,  and  P.M.  Zion,  1987.  Spatial  and  temporal  variability  of  phytoplankton 

pigment  off  Northern  California  during  CODE-1.  J.  Geophy.  Res.  92,  1745-1756. 

Abbott,  M.R.,  and  B.  Barksdale,  1991.  Phytoplankton  Pigment  Patterns  and  Wind  Forcing 
off  Central  California.  J.  Geophys.  Res.,  96,  14,649-14,667. 

Allen,  J.  S.,  L.J.  Walstad,  and  P.A.  Newberger,  1991.  Dynamics  of  the  Coastal  Transition 
Zone  Jet,  2.  Nonlinear  finite  amplitude  behavior,  J.  Geophys.  Res.,  96,  14,995- 
15,016. 

Barth,  J.A.  and  D.M.  Checkley,  2009.  Patterns  and  processes  in  the  California  Current 
System,  Progress  in  Oceanography,  83,  49-64. 

Barth,  J.A.,  S.D.  Pierce  and  R.L.  Smith,  2000.  A  separating  coastal  upwelling  jet  at  Cape 
Blanco,  Oregon  and  its  connection  to  the  California  Current  System.  Deep  Sea 
Research  II,  47,  783-810. 

Belkin,  I.M.,  P.C.  Cornillon  and  K.Sherman,  2009.  Fronts  in  Farge  Marine  Ecosystems. 
Progress  In  Oceanography,  81(1-4),  223-236. 

Bennett,  A.  F.,  1992.  Inverse  Methods  in  Physical  Oceanography.  Cambridge  Un.  Press, 
Cambridge,  346  pp. 

Bennett,  A.  F.,  2002.  Inverse  Modeling  of  the  Ocean  and  Atmosphere.  Cambridge  Un.  Press, 
Cambridge,  234  pp. 

Bernstein,  R.F.,  F.  Breaker,  and  R.  Whitner,  1977.  California  Current  eddy  formation:  Ship, 
air  and  satellite  results.  Science,  195,  353-359. 

Blumberg,  A.F.  and  F.H.  Kantha,  1985.  Open  boundary  condition  for  circulation  models.  J. 
Hydarul.  Eng. -ASCE  111  (2),  237-255. 

Blumberg,  A.F.  and  G.F.  Mellor,  1987.  A  Description  of  a  Three  Dimensional  Coastal  Ocean 
Circulation  Model.  In:  N.S.  Heaps,  Editor,Three-dimensional  coastal  ocean 
models,  American  Geophysical  Union,  Washington  (1987), 1-16. 

Boyer,  T.,  S.  Fevitus,  H.  Garcia,  R.A.  Focarnini,  C.  Stephens  and  J.  Antonov,  2005. 

Objective  analyses  of  annual,  seasonal,  and  monthly  temperature  and  salinity  for 
the  World  Ocean  on  a  0.25°  grid.  International  Journal  of  Climatology,  25(7), 
931-945. 

Brown,  W.  S.,  A.  Gangopadhyay,  F.  F.  Bub,  Z.  Yu  and  G.  Strout,  and  A.  R.  Robinson, 
2007a.  An  Operational  Circulation  Modeling  System  for  the  Gulf  of 
Maine/Georges  Bank  Region,  Part  1:  The  Basic  Elements.  IEEE  J.  Oceanic 
Engineering,  32  (4),  807-822. 

Brown,  W.  S.,  A.  Gangopadhyay,  and  Z.  Yu,  2007b.  An  Operational  Circulation  Modeling 
System  for  the  Gulf  of  Maine/Georges  Bank  Region,  Part  2:  Applications.  IEEE 
J.  Oceanic  Engineering,  32  (4),  807-822. 

Brink,  K.H.,  R.C.  Beardsley,  P.P.  Niiler,  M.  Abbott,  A.  Huyer,  S.  Ramp,  T.  Stanton,  and  D. 

Stuart,  1991.  Statistical  properties  of  near-surface  flow  in  the  California  coastal 
transition  zone.  J.  Geophys.  Res.,  96,  14,693-14,706. 


31 


Brink,  K.H.  and  TJ.  Cowles,  1991.  The  coastal  transition  zone  program.  J.  Geophys.  Res., 

96,  14,637-14,647. 

Brink,  K.  H.,  R.  C.  Beardsley,  J.  Paduan,  R.  Limebumer,  M.  Caruso,  and  J.  G.  Sires,  2000.  A 
view  of  the  1993-1994  California  Current  based  on  surface  drifters,  floats,  and 
remotely  sensed  data.  Journal  of  Geophysical  Research,  105(C4),  8575-8604. 

Bucklin,  A.,  1991.  Population  Genetic  Responses  of  the  Planktonic  Copepod  Metridia 

pacifica  to  a  Coastal  Eddy  in  the  California  Current,  J.  Geophys.  Res.,  96,  14,799- 
14,808. 

Calado,  L,  A.  Gangopadhyay  and  I.C.  daSilveira,  2008:  Feature- Oriented  Regional  Modeling 
and  Simulations  (FORMS)  For  the  Western  South  Atlantic,  Southeastern  Brazil 
Region.  Ocean  Modeling  (25),  48-64. 

Calado,  F.,  A.  Gangopadhyay,  and  I.  C.  A.  da  Silveira,  2006.  A  parametric  model  for  the 
Brazil  Current  meanders  and  eddies  off  southeastern  Brazil,  Geophys.  Res.  Lett, 
33,  F12602,  doi:10.1029/2006GF026092. 

Capet,  X.,  J.  C.  McWilliams,  M.  J.  Molemaker,  A.  F.  Shchepetkin,  2008:  Mesoscale  to 

Submesoscale  Transition  in  the  California  Current  System.  Part  I:  Flow  Structure, 
Eddy  Flux,  and  Observational  Tests.  J.  Phys.  Oceanogr.,  38,  29-43. 

Capet,  X.,  J.  C.  McWilliams,  M.  J.  Molemaker,  A.  F.  Shchepetkin,  2008:  Mesoscale  to 
Submesoscale  Transition  in  the  California  Current  System.  Part  II:  Frontal 
Processes.  J.  Phys.  Oceanogr.,  38,  44-64. 

Castelao,  R.M.,  E.J.D.  Campos  and  J.F.  Miller,  2004.  A  Modelling  Study  of  Coastal 
Upwelling  Driven  by  Wind  and  Meanders  of  the  Brazil  Current.  Journal  of 
Coastal  Research,  (July  2004),  662  -  671. 

Chao,  Y.,  Z.  Fi,  J.  Farrara,  J.  C.  McWilliams,  J.  Bellingham,  X.  Capet,  F.  Chavez,  J.-K.  Choi, 
R.  Davis,  J.  Doyle,  D.  Fratantoni,  P.  Fi,  P.  Marchesiello,  M.  A.  Moline,  J.  Paduan, 
and  S.  Ramp,  2008.  Development,  implementation,  and  evaluation  of  a  data- 
assimilative  ocean  forecasting  system  off  the  central  California  coast.  Deep-Sea 
Research  11,  doi:10.1016/j.dsr2. 2008.08. Oil. 

Chavez,  F.P.  and  C.  Collins,eds,  2000.  Studies  of  the  California  Current  System  Part  2, 
Deep-Sea  Research  11,  47,  5-6. 

Chavez,  F.P.  and  C.A.  Collins,  2000.  Studies  of  the  California  Current  System:  present,  past 
and  future.  Deep-Sea  Research  11,  47,  761-763. 

Chelton,  D.  B.,  1984.  Seasonal  variability  of  alongshore  geostrophic  velocity  off  central 
California.  J.  Geophys.  Res.,  89,  3473-3486. 

Chavez,  F.  P.,  T.J.  Pennington,  R.  Herlien,  H.W.  Jannasch,  G.  Thurmond  and  G.  Friederich, 
1997.  Moorings  and  drifters  for  real-time  interdisciplinary  oceanography,  J. 

Atmos,  and  Ocean.  Tech.,  14,  1199-1211. 

Chavez,  F.P.,  R.T.  Barber,  P.M.  Kosro,  A.  Huyer,  S.R.  Ramp,  T.P.  Stanton,  and  B.  Rojas  de 
Mendiola,  1991.  Horizontal  transport  and  the  distribution  of  nutrients  in  the 
coastal  transition  zone  off  northern  California:  effects  on  primary  production, 
phytoplankton  biomass  and  species  composition.  J.  Geophys.  Res.,  96,  14,833- 
14,848. 


32 


Chelton,  D.  and  M.  Schlax,  1991.  Estimation  of  time  averages  from  irregularly  spaced 
observations:  With  application  to  Coastal  Zone  Color  Scanner  estimates  of 
chlorophyll  concentration.  J.  Geophys.  Res.,  96,  14,669-14,692. 

Chereskin,  T.K.  and  P.P.  Niiler,  1994.  Circulation  in  the  Ensenada  Eront — September  1988. 
Deep  Sea  Research  Part  I:  Oceanographic  Research  Papers,  41,  1251-1287. 

Collins,  C.A.,  J.T.  Pennington,  C.G.  Castro,  T.A.  Rago  and  P.P.  Chavez,  2003.  The 

California  Current  system  off  Monterey,  California:  physical  and  biological 
coupling.  Deep-Sea  Research  11  50  (2003)  2389-2404 

Collins,  C.A.,  C.G.  Castro,  H.  Asanuma,  T.A.  Rago,  H.-K.  Han,  R.  Durazo,  and  P.P.  Chavez, 
2002.  Changes  in  the  hydrography  of  Central  California  waters  associated  with 
the  1997-98  El  Nino.  Progr.  in  Oceanog.,  54,  184-204. 

Collins,  C.A.,  N.  Garfield,  T.A.  Rago,  P.W.  Rischmiller,  and  E.  Carter,  2000.  Mean  Structure 
of  the  Inshore  Countercurrent  and  California  Undercurrent  off  Point  Sur, 
California.  Deep  Sea  Res.,  47,  765-782. 

Collins,  C.  A.,  N.  Garfield,  R.G.  Paquette,  and  E.  Carter,  1996.  Eagrangian  Measurement  of 
subsurface  poleward  Plow  between  380N  and  430N  along  the  West  Coast  of  the 
United  States  during  Summer,  1993.  Geophys.  Res.  Letters,  September  1,  1996, 
23,  No.  18,2461-2464. 

Cornuelle,  B.D.,  T.K.  Chereskin,  P.P.  Niiler,  M.Y.  Morris,  and  D.E.  Musgrave,  2000. 

Observations  and  Modeling  of  a  CA,  Undercurrent  Eddy,”  J.  Geophys.  Res.,  C 
105,  1227-1243. 

Crawford  G.  and  S.  Stone,  2008.  COCMP  Surface  Current  Mapping  Reveals.  Eddy  and 

Upwelling  Jet  off  Cape  Mendocino.  American  Geophysical  Union,  Pall  Meeting 
2008,  abstract  #OS41B-1216. 

Cummings,  J.A.,  C.  Szczechowski  and  M.  Carnes,  1997.  Global  and  regional  ocean  thermal 
analysis  systems.  Marine  Tech.  Soc.  J.,  31,  63-75. 

Dewey,  R.K.,  J.N.  Mourn,  C.A.  Paulson,  D.R.  Caldwell,  and  S.D.  Pierce,  1991.  Structure  and 
dynamics  of  a  coastal  filament,  J.  Geophys.  Res.,  96,  14,885-14,907. 

Ezer,  T  and  G.E.  Mellor,  1994.  Diagnostic  ands  prognostic  calculations  of  the  North  Atlantic 
circulation  and  sea  level  using  a  sigma  coordinate  ocean  model.  J.  Geophy.  Res., 
99,  14,159-  14,171. 

Ezer  T.  and  G.E.  Mellor,  1997.  Data  assimilation  experiments  in  the  Gulf  Stream  region: 

how  useful  are  satellite-derived  surface  data  for  nowcasting  the  subsurface  fields? 
J  Atmos  Ocean  Tech,  16(6),  1379-1391. 

Plament,  P.,  E.  Armi,  and  E.  Washburn,  1985.  The  evolving  structure  of  an  upwelling 
filament.  J.  Geophys.  Res.,  90:1765-1985. 

Pox  D.N.,  M.R.  Carnes  and  J.E.  Mitchell,  1992.  Characterizing  Major  Prontal  Systems:  A 
Nowcast  /Porecast  System  for  Northwest  Atlantic.  Oceanography,  5,  49-53. 

Pratantoni,  D.M.  and  S.  H.  D.  Haddock,  2009.  Introduction  to  the  Autonomous  Ocean 

Sampling  Network  (AOSN)  program.  Deep-Sea  Research  II,  Deep-Sea  Research 
77,56,61. 


33 


Gangopadhyay,  A.,  A.R.  Robinson,  and  H.G.  Arango,  1997.  Circulation  and  Dynamics  of  the 
Western  North  Atlantic,  I:  Multi-scale  Feature  Models.  J.  Atmos,  and  Oceanic 
Tech.,  14(6),  1314-1332 

Gangopadhyay,  A.  and  A.R.  Robinson,  2002.  Feature- Oriented  Regional  Modeling  of 
Oceanic  Fronts.  Dynamics  of  Atmospheres  and  Oceans,  36,  201-232. 

Gangopadhyay,  A.,  A.R.  Robinson,  P.J.  Haley,  Jr.,  W.G.  Leslie,  C.J.  Lozano,  J.J.  Bisagni  and 
Z.  Yu,  2003.  Feature  Oriented  Regional  Modeling  and  Simulations  (FORMS)  in 
the  Gulf  of  Maine  and  Georges  Bank.  2003,  Continental  Shelf  Research.  23  (3-4), 
317-353. 

Garfield,  N.,  M.E.  Maltrud,  C.A.  Collins,  T.A.  Rago  and  R.G.  Paquette,  2001.  Lagrangian 
flow  in  the  California  Undercurrent,  an  observation  and  model  comparison. 
Journal  of  Marine  Systems,  29(1-4),  201-220. 

Garfield,  N.,  C.A.  Collins,  R.G.  Paquette,  and  E.  Carter,  1999.  Eagrangian  Exploration  of  the 
California  Undercurrent,  1992-1995.  J.  Phys.  Oceanogr.,  29,  560-583. 

Glenn,  S.M.,  and  A.R.  Robinson,  1995.  Validation  of  an  operational  Gulf  Stream  forecasting 
model.  Qualitative  Skill  assessment  for  coastal  models,  AGU  Estuarine/Coastal 
Series,  Vol.  47,  American  Geophysical  Union,  469-499. 

Graf,  J.,  Sasaki,  C.,  Winn,  C.,  Eiu,  W.  T.,  Tsai,  W.,  Ereilich,  M.,  Eong,  D.,  1998.  NASA 
scatterometer  experiment.  Acta  Astronaut.  43,  397-407. 

Haidvogel,  D.  U.,  A.  Beckmann,  and  K.  S.  Hedstrom.  1991.  Dynamical  simulations  of 

filament  formation  and  evolution  in  the  Coastal  Transition  Zone.  J.  Geophys.  Res. 
96:15,017-15,040. 

Haley  Jr.,  P.J.,  P.E.J.  Ixrmusiaux,  A.R.  Robinson,  W.  Eeslie,  O.  Eogoutov,  G.  Cossarini, 

X.S.  Eiang,  P.  Moreno,  S.R.  Ramp,  J.  Doyle,  J.G.  Bellingham,  E.P.  Chavez,  and 
S.  Johnston  (2009).  Porecasting  and  reanalysis  in  the  Monterey  Bay/California 
current  region  for  the  Autonomous  Ocean  Sampling  Network- II  experiment. 
Deep-Sea  Research  Part  II,  56:  127-148. 

Haynes,  R.,  and  E.D.  Barton,  1991.  Eagrangian  observations  in  the  Iberian  coastal  transition 
zone,  J.  Geophys.  Res.,  96,  14,731-14,742. 

Hayward,  T.E.  and  A.W.  Mantyla,  1990.  Physical,  chemical  and  biological  structure  of  a 
coastal  eddy  near  Cape  Mendocino,  J.  Mar.  Res.,  4,  825-850. 

Hickey,  B.M.,  1998.  Coastal  Oceanography  of  Western  North  America  from  the  tip  of  Baja 
California  to  Vancouver  Is.,  pp  345-393  in  Volume  11,  Chapter  12,  THE  SEA, 
eds.  K.H.  Brink  and  A.R.  Robinson,  Wiley  and  Sons,  Inc. 

Hulburt  H.  E.,  D.N.  Pox  and  E.J.  Metzger,  1990.  Statistical  Inference  of  Weakly  Correlated 
Subthermocline  Pields  Prom  Satellite  Altimeter  Data.  J.  Geophy.  Res.,  95, 
11,375-11,409. 

Hurlburt,  H.E.,  A.J.  Wallcraft,  W.J.  Schmitz,  Jr.,  P.J.  Hogan  and  E.J.  Metzger,  1996. 

Dynamics  of  the  Kuroshio/Oyashio  current  system  using  eddy-resolving  models 
of  the  North  Pacific  Ocean,  J.Geophys.Res.,  101,  941-976. 


34 


Huyer,  A.E.,  1983.  Coastal  upwelling  in  the  California  Current  system.  Progr.  Oceanogr.  12, 
259-284. 

Huyer,  A.,  P.M.  Kosro,  J.  Fleischbein,  S.R.  Ramp,  T.  Stanton,  L.  Washburn,  F.P.  Chavez, 

T.J.  Cowles,  S.D.  Pierce,  and  R.F.  Smith,  1991.  Currents  and  water  masses  of  the 
Coastal  Transition  Zone  off  northern  California,  June  to  August,  1988,  J. 

Geophys.  Res.,  96,  14,809-14,831. 

Huyer,  A.,  J.  A.  Barth,  P.  M.  Kosro,  R.  K.  Shearman,  and  R.  F.  Smith,  1998.  Upper-ocean 
water  mass  characteristics  of  the  California  Current,  summer  1993,  Deep  Sea 
Res.,  45,  Part  II,  1411-1442. 

Ikeda,  M.  and  W.J.  Emery,  1984.  Satellite  observations  and  modeling  of  meanders  in  the 
California  Current  System  off  Oregon  and  Northern  California,  J.  Phys. 
Oceanography,  14,  1434-1450. 

Kadko,  D.C.,  F.  Washburn,  and  B.  Jones.  1991.  Evidence  of  subduction  within  cold 

filaments  of  the  northern  California  coastal  transition  zone.  J.  Geophys.  Res.,  96, 
14,909-14,926. 

Kelly,  K.  A.,  R.  C.  Beardsley,  R.  Fimeburner,  K.  H.  Brink,  J.  D.  Paduan,  and  T.  K.  Chereskin 
(1998),  Variability  of  the  near-surface  eddy  kinetic  energy  in  the  California 
Current  based  on  altimetric,  drifter,  and  moored  current  data,  J.  Geophys.  Res., 
103(C6),  13,067-13,083,  doi:10.1029/97JC03760. 

Kim,  H-S.,  A.  Gangopadhyay,  F.  K.  Rosenfeld,  and  F.F.  Bub,  2007:  A  High-resolution 
Regional  Climatology  for  the  Central  California,  Continental  Shelf  Res. 
doi:10.1016/j.csr.2007.05.011. 

Kosro,  P.M.,  A.  Huyer,  S.R.  Ramp,  R.F.  Smith,  F.P.  Chavez,  T.F.  Cowles,  M.R.  Abbott,  P.T. 

Strub,  R.T.  Barber,  P.  lessen  and  F.F.  Small,  1991.  The  structure  of  the  transition 
zone  between  coastal  waters  and  the  open  ocean  off  northern  California,  winter 
and  spring  1987.  J.  Geophys.  Res.,  96,  14,707-14,730. 

Fagerloef,  G.S.E.,  1992.  The  Point  Arena  Eddy:  a  recurring  summer  anticyclone  in  the 
California  Current,  J.  Geophys.  Res.,  97,  C8,  12,557-12,568. 

Fermusiaux,  P.F.J.,  1999a.  Data  Assimilation  via  Error  Subspace  Statistical  Estimation,  Part 
II:  Mid- Atlantic  Bight  Shelfbreak  Front  Simulations  and  ESSE  Validation.  Mon 
Weath.  Rev,  127(8),  1408-1432. 

Fermusiaux,  P.,  1999b.  Estimation  and  study  of  mesoscale  variability  in  the  strait  of  Sicily. 
Dynam.  Atmos.  Oceans  29,  255-303. 

Fozano,  C.J.,  A.R.  Robinson,  H.G.  Arango,  A.  Gangopadhyay,  N.Q.  Sloan,  P.J.  Haley,  F. 

Anderson,  and  W.G.  Feslie,  1996.  An  Interdisciplinary  Ocean  Prediction  System: 
Assimilation  Strategies  and  Structured  Data  Models.  In:  Modem  Approaches  to 
Data  Assimilation  on  Ocean  Modeling,  P.  Mallanotte-Rizzoli  (ed.),  Elsevier 
Oceanography  Series,  61,  Amsterdam,  pp.  413-452. 

Fynn,  R.J.,  K.A.  Bliss  and  F.F.  Eber,  1982.  Vertical  and  horizontal  distributions  of  seasonal 
mean  temperature,  salinity,  sigma-t,  stability,  dynamic  height,  oxygen,  and 
oxygen  saturation  in  the  California  Current,  1950-1978.  CalCOFI  Atlas  No.  30, 
Scripps  Institution  of  Oceanography,  Fa  Jolla,  513 


35 


Lynn,  R.  J.,  and  JJ.  Simpson,  1987.  The  California  Current  system:  the  seasonal  variability 
of  its  physical  characteristics,  J.  Geophys.  Res.,  92,  12,947-12,966. 

Mackas,  D.L.,  L.  Washburn,  and  S.L.  Smith.  1991.  Zooplankton  community  pattern 

associated  with  a  California  Current  cold  filament.  J.  Geophys.  Res.,  96,  14,781- 
14,797. 

Magnell,  B.A.,  N.A.  Bray,  C.D.  Winant,  C.L.  Greengrove,  J.  Largier,  J.F.  Borchard,  R.L. 
Bernstein,  C.E.  Dorman,  1990.  Convergent  shelf  flow  at  Cape  Mendocino. 
Oceanography,  3(1),  4-11. 

Marchesiello,  P.,  J.C.  McWilliams,  and  A.  Shchepetkin,  2003.  Equilibrium  structure  and 
dynamics  of  the  California  Current  System.  J.  Phys.  Oceanogr.,  33,  753-783. 

Mathew  E.M.  and  J.E.  McClean,  2005.  An  eddy  resolving  global  1/10°  ocean  simulation. 
Ocean  Modelling,  8,  31  -  54. 

McWilliams,  J.C.,  1985.  Submesoscale,  coherent  vortices  in  the  ocean.  Rev.  of  Geophys.,  23, 
165-182. 

Mied,  R.P.,  1990.  Mushroom-Eike  Patterns  of  the  Ocean  Surface,  Eos,  Transactions, 
American  Geophysical  Union,  71  (49),  1837. 

Mied,  R.P.,  J.C.  McWilliams,  and  G.J.  Eindemann,  1991.  The  generation  and  evolution  of 
mushroom-like  vortices,  J.  Phys.  Ocean.,  21,  490-510. 

Miller,  A.J.,  J.C.  McWilliams,  N.  Schneider,  J.S.  Allen,  J.A.  Barth,  R.C.  Beardsley,  E.P. 
Chavez,  T.K.  Chereskin,  C.A.  Edwards,  R.E.  Haney,  K.A.  Kelly,  J.C. 

Kindle,  E.N.  Ey,  J.R.  Moisan,  M.A.  Noble,  P.P.  Niiler,  E.Y.  Oey,  E.B.  Schwing, 
R.K.  Shearman,  and  M.S.  Swenson.  1999.  Observing  and  Modeling  the  California 
Current  System,  Eos,  Transactions,  American  Geophysical  Union,  80,  533-539. 

Miller,  A.J.,  E.  Di  Eorenzo,  D.J.  Neilson,  B.D.  Cornuelle  and  J.R.  Moisan,  2000.  “Modeling 
CalCOEI  Observations  during  El  Nino:  Fitting  physics  and  biology,”  California 
Cooperative  Oceanic  Fisheries  Investigations  Reports,  41,  87-97. 

Moore,  A.  M.,  Hernan  G.  Arango,  Emanuele  Di  Eorenzo,  Arthur  J.  Miller,  Bruce  D. 

Cornuelle,  2009:  An  Adjoint  Sensitivity  Analysis  of  the  Southern  California 
Current  Circulation  and  Ecosystem.  J.  Phys.  Oceanogr.,  39,  702-720. 

Mooers,  C.N.K.  and  A.R.  Robinson,  1984.  Turbulent  Jets  and  Eddies  in  the  California 
Current  and  Inferred  Cross-Shore  Transports.  Science,  223,  51-53. 

Munk,  W.,  E.  Armi,  K.  Eischer.  and  E.  Zachariasen,  2000.  Spirals  on  the  sea,  Proc.  R.  Soc. 
Lond.,  456,  1217-1280. 

Noble,  M.  A.  and  S.  R.  Ramp,  2000.  Subtidal  currents  over  the  central  California  slope: 

Evidence  for  offshore  veering  of  the  undercurrent  and  for  direct,  wind-driven 
slope  currents.  Deep-Sea  Res.  II,  47,  871-906. 

Oey,  E.,  1999.  A  forcing  mechanism  for  the  poleward  flow  off  the  southern  California  coast, 
J.  Geophys.  Res.,  104,  13,529-13,539. 

Palma,  E.  D.,  Matano,  R.  P.,  1998.  On  the  implementation  of  passive  open  boundary 

conditions  for  a  general  circulation  model:  The  barotropic  mode.  J.  Geophys.  Res. 
103,  1319-1341. 


36 


Penven,  P.,  Laurent  Debre,  Patrick  Marchesiello  and  J.  C.  McWilliams,  2006.  Evaluation  and 
application  of  the  ROMS  1-way  embedding  procedure  to  the  central  California 
upwelling  svstem,Ocean  Modelling,  Volume  12,  Issues  1-2,  2006,  Pages  157-187. 

Pierce,  S.  D.,  R.L.  Smith,  P.M.  Kosro,  J.A.  Barth,  and  C.D.  Wilson,  2000.  Continuity  of  the 
poleward  undercurrent  along  the  eastern  boundary  of  the  mid-latitude  Pacific, 
Deep  Sea  Res.,  47,  811-829. 

Pringle,  J.M.,  Edward  P.  Dever.  (2009)  Dynamics  of  wind-driven  upwelling  and  relaxation 
between  Monterey  Bay  and  Point  Arena:  Eocal-,  regional-,  and  gyre-scale 
controls.  Journal  of  Geophysical  Research  114:C7, 

Ramp,  S.R.,  J.E.  McClean,  C.A.  Collins,  A.J.  Semtner  and  K.A.S.  Hays,  1997a.  Observations 
and  modeling  of  the  1991-1992  El  Nino  signal  off  central  California,  J.  Geophys. 
Res.,  102,  5553-5582. 

Ramp,  S.R.,  E.K.  Rosenfeld,  T.D.  Tisch,  and  M.R.  Hicks,  1997b.  Moored  observations  of 
the  current  and  temperature  structure  over  the  continental  slope  off  central 
California,  1,  A  basic  description  of  the  variability,  J.  Geophys.  Res.,  102, 
22,877-22,902. 

Ramp,  S.R.,  P.E.  lessen,  K.H.  Brink,  P.P.  Niiler,  E.E.  Daggett  and  J.S.  Best,  1991.  The 
physical  structure  of  cold  filaments  near  Point  Arena,  California,  During  June 
1987.  J.  Geophys.  Res.,  96„  14,859-14,883. 

Ramp  S.R.,  J.D.  Paduan,  I.  Shulman,  J.  Kindle,  P.E.  Bahr  and  P.  Chavez,  2005.  Observations 
of  upwelling  and  relaxation  events  in  the  northern  Monterey  Bay  during  August 
2000.  J.  Geophys.  Re5.,110,  C07013,  doi:10.1029/2004JC002538. 

Ramp,  S.  R.,  R.  E.  Davis,  N.  E.  Eeonard,  I.  Shulman,  Y.  Chao,  A.  R.  Robinson,  J.  Marsden, 

P.  P.  J.  Eermusiaux,  D.  Pratantoni,  J.  D.  Paduan,  P.  P.  Chavez,  P.  E.  Bahr,  S. 
Eiang,  W.  Eeslie,  and  Z.  Pi,  2008.  Preparing  to  predict:  The  Second  Autonomous 
Ocean  Sampling  Network  (AOSN-II)  experiment  in  the  Monterey  Bay.  Deep-Sea 
Research  II,  doi:  10. 1016/j.dsr2. 2008.08. 013. 

Rischmiller,  P.W.,  1993.  Variability  of  the  California  Current  system  off  Point  Sur, 
California,  from  April  1988  to  December  1990.  M.S.  Thesis,  Naval  Postgraduate 
School,  Monterey,  157. 

Robinson,  A.  R.,  S.M.  Glenn,  M.A.  Spall,  E.J.  Walstad,  G.M.  Gardner,  and  W.G. Eeslie, 

1989.  Porecasting  meanders  and  Rings,  EOS.  Oceanogr.  Report,  70(45),  1464- 
1473. 

Robinson,  A.  R.  and  S.M.  Glenn,  1999.  Adaptive  Sampling  for  Ocean  Porecasting.  Naval 
Research  Reviews,  51(2),  28-38. 

Robinson,  A.R.,  J.A.  Carton,  C.N.K.  Mooers,  E.J.  Walstad,  E.E.  Carter,  M.M.  Rienecker, 
J.A.  Smith,  and  W.G.  Eeslie,  1984.  A  Real-Time  Dynamical  Porecast  of  Ocean 
Synoptic/Mesoscale  Eddies.  Nature,  309(5971),  781-783. 

Robinson,  A.R.,  J.A.  Carton,  N.  Pinardi  and  C.N.K.  Mooers,  1986.  Dynamical  Porecasting 
and  Dynamical  Interpolation: An  Experiment  in  the  California  Current.  Journal  of 
Physical  Oceanography,  16,  1561-1579. 


37 


Robinson,  A.R.,  and  A.  Gangopadhyay,  1997.  Circulation  and  Dynamics  of  the  Western 
North  Atlantic,  II:  Dynamics  of  Meanders  and  Rings.  Journal  of  Atmospheric  and 
Oceanic  Technology,  14(6),  1333-1351. 

Rosenfeld,  L.K.,  F.B.  Schwing,  N.  Garfield,  and  D.E.  Tracy,  1994.  Bifurcated  flow  from  an 
upwelling  center:  a  cold  source  for  Monterey  Bay.  Continental  Shelf 
Research,  14,  931-964. 

Schofield,  O.,  S.M.  Glenn,  J.  Orcutt,  M.  Arrott,  M.  Meisinger,  A.  Gangopadhyay,  W.  Brown, 
R.  Signell,  M.  Moline,  Y.  Chao,  S.  Chien,  D.  Thompson,  A.  Balasuriya,  P. 
Lermusiaux  and  M.  Oliver,  2010:  Automated  Sensor  Networks  to  Advance  Ocean 
Science,  EOS,  91  (39),  28  Sept,  2010. 

Shaji,  C.  and  A.  Gangopadhyay,  2007.  Synoptic  modeling  of  the  West  India  Coastal  Current 
System  using  an  upwelling  feature  model.  Journal  of  Atmospheric  and  Oceanic 
Technology,  24(5),  877-893. 

Sheres,  D.  and  K.E.  Kenyon,  1989.  A  double  vortex  along  the  California  Coast,  J.  Geo.  Res., 
94,  4989-4997. 

Shulman,  L,  C.  Rowley,  S.  Anderson,  S.  DeRada,  J.  Kindle,  P.  Martin,  J.  Doyle,  J. 

Cummings,  S.  Ramp,  E.  Chavez,  D.  Eratantoni,  and  R.  Davis,  2008.  Impact  of 
glider  data  assimilation  on  the  Monterey  Bay  model.  Deep-Sea  Research  II, 
doi:10.1016/j.dsr2.2008.08.003. 

Smith,  S.E.,  and  P.V.Z.  Eane,  1991.  The  jet  off  Point  Arena,  California:  its  role  in  aspects  of 
secondary  production  in  the  copepod  Eucalanus  californicus  Johnson.  J.  Geophys. 
Res.,  96,  14,849-14,858. 

Strub,  P.T.,  P.M.  Kosro  and  A.  Huyer,  1991.  The  nature  of  the  cold  filaments  in  the 
California  Current  system.  J.  Geophys.  Res.,  96,  14,743-14,768. 

Strub,  P.  T.  and  C.  James,  2000.  Altimeter-derived  variability  of  surface  velocities  in  the 

California  Current  System:  2.  Seasonal  circulation  and  eddy  statistics.  Deep-Sea 
Res.  II,  45,  831-870. 

Swenson,  M.  S.,  and  P.  P.  Niiler,  1996.  Statistical  analysis  of  the  surface  circulation  of  the 
California  Current,  J.  Geophys.  Res.,  101,  22,631-22,645. 

Tisch,  T.  D.,  S.R.  Ramp,  and  C.A.  Collins,  1992.  Observations  of  the  geostrophic  current  and 
water  mass  characteristics  off  Point  Sur,  California,  from  May  1988  through 
November  1989.  J.  Geophys.  Res.,  97,  12,535-12,555. 

Tracy,  D.,  1990.  Source  of  cold  water  in  Monterey  Bay  observed  by  AVHRR  satellite 
imagery,  MS  thesis.  Naval  Postgraduate  School,  Monterey,  CA. 

Traganza,  E.D.,  D.A.  Nestor,  and  A.K.  McDonald,  1980.  Satellite  observations  of  a  nutrient 
upwelling  off  the  coast  of  California,  J.  Geophys.  Res.,  85,  4,101-4,106. 

Traganza,  E.D.,  J.C.  Conrad,  and  E.C.  Breaker,  1981.  Satellite  observations  of  a  cyclonic 
upwelling  system  and  giant  plume  in  the  California  Current,  pp.  229-241  in 
Coastal  Upwelling,  E.  A.  Richards  (ed),  American  Geophysical  Union, 
Washington,  DC. 


38 


Walstad,  L.J.,  J.S.  Allen,  P.M.  Kosro,  and  A.  Huyer,  1991.  Dynamics  of  the  Coastal 

Transition  Zone  through  data  assimilation  studies,  J.  Geophys.  Res.,  96,  14,959- 
14,978. 

Warn-Varnas,  A.,  A.  Gangopadhyay,  and  J.  Hawkins,  2007:  Water  Masses  in  the  Monterey 
Bay  during  summer  of  2000.  Continental  Shelf  Research, 
doi:10.1016/j.csr.2007.01.004. 

Washburn,  L.,  D.C.  Kadco,  B.H.  Jones,  T.  Hayward,  P.M.  Kosro,  T.P.  Stanton,  S.  Ramp,  and 
T.  Cowles,  1991.  Water  mass  subduction  and  the  transport  of  phytoplankton  in  a 
coastal  upwelling  system.  J.  Geophys.  Res.,  96,  14,927-14,945. 

Wickham,  J.B.,  1975.  Observations  of  the  California  countercurrent.  J.  Mar.  Res.,  33,  325- 
340. 

Wickham,  J.  B.,  A.A.  Bird,  and  C.N.K.  Mooers,  1987.  Mean  and  variable  flow  over  the 
central  California  continental  margin.  Cont.  Shelf  Res.,  7,  827-849. 

Wooster,  W.S.  and  J.H.  Jones,  1970.  California  Undercurrent  off  Northern  Baja  California,  J. 
Mar.  Res.,  28,  235-250. 

Wyllie,  J.  G.  1966.  Geostrophic  flow  of  the  California  Current  at  the  surface  and  at  200 
meters.  Calif.  Coop.  Oceanic  Fish.  Invest.,  Atlas  No.4,  288  pp. 


39 


Table 


Table  1.  Characteristics  of  the  synoptic  features  for  the  California  Current  System. 
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Location 

Width 
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Vertical 
extent  (m) 
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Variability 

Represented 

by 

California 

Current 

(southward 

flow) 

100-150 
km  offshore 
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1350 

0-500 

Baroclinic  jet, 

Vmax  =  50-70 
cm/s.  Wavelength 
=  0  (300km)  Time- 
scale  =  5-10  days 

Mean  flow 
(wide  front) 
or 

climatology 

CUC/Inshore 

Current 

(poleward 

flow) 

10  -  100  km 
offshore 

10-40 

0-300 

Core  @  100-200  m 
depths: 

Vmax  = 

15-20  cm/s 

Low-salinity 

Parabolic 

pool 

Coastal 

Transition 

Zone 

Continental 

shelf 

Filaments,  eddies, 
upwelling  fronts 
and  anomalous 
pools 

Water- mass 
melding 

Coastal 

eddies 

C. 

Mendocino 
and  Pt. 

Arena 

<  100 
km 

<300 

Low  T,  hi  S  and  hi 
Nutrient; 

Time-scale  =  40- 
60  days 

Asymmetric 

eddies 

Coastal  jets 
along 
upwelling 
fronts 

Over  the 
shelf 

10- 

40km 

Above 

halocline 

(00=26.5) 

Vmax  =  100  cm/s, 
O(week) 

Frontal 

water-mass 

structure 

Anomalous 

pools 

Onshore 
side  of 
upwelling 
fronts 

20- 

30km 

<60 

Low  T  and  low  S, 
Time-scale  = 
C(week) 

Eddy 

structure 

Large 

filaments 

Up  to  200 
km  offshore, 
and  near  Pt. 
Arena 

< 

100km 

100  or 
below 

LowT(12-13C) 
and  hi  S  (32.7- 
33psu),  Time-scale 
=  2-4  weeks; 
Vertical  velocity  = 
-25  to  40  m/day 

Weak  frontal 

structure 

Smaller 

niaments 

(squirts) 

Onshore  of 

upwelling 

fronts 

30  and 
50-100 
km 
long 

Near 

surface 

LowT(10-12C) 
and  hi  S  (>  33psu); 
Time-scale  =  6-10 
days 

Strong  frontal 
signature 

Mushroom 

heads 

T- shape 

Above 

seasonal 

thermocline 

Ageostrophic  & 
asymmetric. 
Time-scale  =  1-5 
days 

Asymmetric 
eddies  and  a 
front 

Features 

References 

California  Current 
(Mean  and  its 
baroclinic  jet  core) 

Brink  et  ah,  1991;Lynn  and  Simpson,  1987;  Chelton,  1984;  Tisch  et  al., 
1992;  Chereskin  et  al.,  1998; 

Miller  et  al.,  1999;  Ramp  et  al.,  1997a,b;  Collins  et  al.,  2003 

Poleward  Flow 
(California 
Undercurrent  and 
Inshore  Current 
(Davidson 

Current)) 

Ramp  et  al.,  1997a,b;  Wooster  and  Jones,  1970; 

Wickham,  1975;  Pierce  et  al.,  2000;  Oey,  1999;  Swenson  and  Niiler, 
1996;  Huyer  et  al.,  1998;  Chavez  et  al.,  1997,  Garfield  et  al.,  1999; 
Collins  et  al.,  1996;  Marchesiello  et  al.,  2003;  Collins  et  al.,  2000,  2002, 

2003 

Coastal  Transition 
Zone  (CTZ) 

Brink  and  Cowles,  1991;  Brink  et  al.,  1991;  Chelton  and  Schlax,  1991; 
Kosro  et  al.,  1991;  Haynes  and  Barton,  1991;  Strub  et  al.,  1991; 
Marchesiello  et  al.,  2003 

Coastal  eddies 

Hayward  and  Mantyla,  1990;  Bucklin,  1991;  Hickey,  1998. 

Coastal  jets  along 
Upwelling  fronts 

Huyer  et  al.,  1991;  Smith  and  Lane,  1991;  Pierce  et  al.,  1991;  Allen  et 
al.,  1991;  Kosro  et  al.,  1991;  Strub  et  al.,  1991;  Rosenfeld  et  a/.,  1994; 
Chavez  et  al.,  1997;  Huyer,  1983;  Washburn  et  al.,  1991 

Anomalous  pools 

Hayward  and  Mantyla,  1990;  Strub  et  al.,  1991 

Large  filaments 

Bernstein  et  al.,  1977;  Traganza  et  al.,  1980,  1981;  Flament  et  al.,  1985; 
Abbott  and  Zion,  1987;  Brink,  1991;  Strub  et  al.,  1991;  Mackas  etal., 
1991;  Ramp  et  a/.,  1991;  Dewey  etal.,  1991;  Kadko  eta/.,  1991; 
Chavez  et  al.,  1991;  Chereskin  and  Niiler,  1994 

Squirts  (smaller 
filaments) 

All  of  the  above,  specially  Ramp  et  al.,  1991;  Dewey  et  al.,  1991; 

Hickey,  1998 

Mushroom  heads 

Ikeda  and  Emery,  1984;  Sheres  and  Kenyon,  1989;  Smith  et  al.,  1991; 
Mied,  1990;  Mied  et  al.,  1991;  Munk  et  al.,  2000 

Table  2.  Reference  list  for  the  studies  that  were  used  to  select  parameter  ranges  for  the  features. 


Saiinity 

(psu) 

Max  core 
width(km) 

Depth(km) 

a(km) 

Line  60 

33 

180.92105 

0.09 

22730.85 

33.1 

340.13158 

0.103 

70199.93 

33.2 

455.92105 

0.115 

112969.6 

33.3 

557.23684 

0.127 

152811.5 

33.4 

651.31579 

0.139 

190742.9 

33.5 

745.39474 

0.151 

229972.4 

33.6 

832.23684 

0.163 

265574.4 

33.7 

940.78947 

0.175 

316101.7 

Line  90 

33.3 

79.605263 

0.068 

5824.447 

33.4 

405.26316 

0.1 

102648.9 

33.5 

506.57895 

0.126 

127292.8 

33.6 

607.89474 

0.15 

153973.3 

33.7 

723.68421 

0.174 

188117.4 

Line  110 

33.6 

383.55263 

0.11 

83586.72 

33.7 

528.28947 

0.12 

145359.3 

33.8 

644.07895 

0.135 

192054.5 

Table  3.  Parameter  values  for  the  low-salinity  pool  parabolic  feature  model. 


Table  4.  Nature  and  Characteristics  of  Simulated  vs.  observed  eddies  in  the 
CCS  region  (36-48N) 


Eddy  Type 

Eddy 

Center 

Size 

(simulated) 

Size 

(Observed) 

Observational 

Evidence 

(Reference) 

Comments 

ACEl 

Pt  Sur 

25-50  km 

R  of  50- 
100  km 

Ramp  et  al. 
(1997b) 
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flow 
reversal 
down  to 
1000m 
(>20cm/sec 
in  300m; 
>30cm  at 
100m) 

ACE2 

Pt  Arena 

50  km 

R  of  50- 
100  km 

Eagerloef 

(1992) 

CW-CCW 

pairs 

Off 

Monterey 

Bay 

>50  km 

50-100  km 

Eagerloef 

1992  (Eig.  9) 

Cuddies 

cue 

Eddy 

25-35  km 

Comuelle  et 
al.  2000 

Anticyclones 

Cuddies 

At 

different 
locations 
near  Capes 

R  <  35  km 

Garfield  et  al. 
1999,  2001 

Propagate 
westward 
at  1-2 
cm/sec 

Earger 
eddies  at 
300m 

R>25km 

R>  50km 

Garfield  et  al. 
1999,  2001 

Figure 
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Figure  1(a).  A  typical  SST  for  the  California  Current  System  shows  the  complex  multiscale 
phenomena  along  this  eastern  boundary  current  region,  (b)  A  typical  climatological  surface 
temperature  field  for  August  (from  Levitus). 
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Figure  2. 

The  three  simple  conceptual  models  considered  for  the  flow  structure  associated  with  the  cold 
filaments  as  proposed  by  Strub  et  al.  (1991):  Squirts  (left),  mesoscale  eddy  field  (middle)  and 
meandering  jet  (right). 
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Figure  3. 

This  schematic  shows  the  complexity  in  the  meander-eddy-upwelling  system  of  the  California 
Current  off  Oregon  and  California  in  the  spring/summer.  The  coastal  jet  coming  from  the  north 
moves  offshore  as  the  season  progresses  (denoted  by  the  March,  May,  and  July  lines),  and 
contributes  downstream  to  the  low-salinity  core  of  the  California  Current.  This  core  delineates 
the  higher  variability  region  nearshore  from  the  less  active  regions  offshore.  The  offshore  region 
consists  of  the  mean  southeastward  flow  of  the  wind-driven  subtropical  gyre.  The  inshore  region 
is  populated  by  upwelling  centers  concentrated  near  the  capes.  Eddies  generated  by  the 
meandering  jet  are  found  in  both  regions,  and  may  exhibit  either  cyclonic  or  anticyclonic 
rotation.  Temporally  transient  eddies  of  various  sizes  have  been  observed. 
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Figure  4.  (a)  The  CalCOFI  station  locations  along  different  lines  [from  Lynn  and  Simpson, 
1987];  (b)  (c)  and  (d)  are  the  vertical  sections  of  mean  salinity  for  the  month  of  July  along 
CalCOFI  station  lines  60,  90  and  1 10  [from  Lynn  et  ah,  1982]  Note  the  consistency  of  the  low 
salinity  pool  in  these  sections. 
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Figure  5.  (a)  Schematic  representation  of  a  parabolic  isohaline  for  the  low-salinity  pool  feature 
model;  (b)  Feature  modeled  vertical  section  of  salinity  at  line  60. 
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Figure  6.  The  temperature  and  salinity  profiles  for  the  California  Current  Core. 
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Figure.  7.  Locations  of  ADCP  sections  and  maximum  speed  and  locations  of 
the  California  Undercurrent.  (From  Pierce  et  al.  2000).  ADCP  sections  were 
made  across  the  continental  margin  from  32-5  IN  at  about  20  km  spacing  during 
the  1995  upwelling  season.  The  depth  of  the  core  of  poleward  flow  at  each 
section  is  determined  by  taking  the  center  of  mass  (these  depths  range  between 
118  and  258  m).  The  vertical  arrows  show  the  speed  and  location  of  the 
maximum  poleward  flow  found  at  this  depth.  The  thick  horizontal  lines  show 
the  width  of  the  poleward  flow  at  half-maximum  velocity.  The  whiskers  show 
the  full  lengths  of  the  ADCP  sections. 


7 


Temperature  Salinity 


Figure  8.  The  CUC-FM  (black)  profiles  compared  with  observed  (grey)  and  Levitus  (color) 
profiles. 
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Figure  9.  A  schematic  of  the  upwelling  feature  model  (after  Calado  et  ah,  2008).  See  text  for  the 
details  of  the  parameter  definition. 
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CUF  FM  T-profile  @Cape  Blanco 


CUF  FM  S-profile  @Cape  Bianco 


Figure  10.  Coastal  upwelling  front  feature  model  T-profiles  (left  —  a)  and  S-profiles  (right  —  b) 
at  Cape  Blanco.  Thick  blue  and  red  lines  represent  inshore  and  offshore  profiles,  superimposed 
on  generated  FM  T&S  along  a  transect.  The  bottom  panel  (c)  shows  the  T/S  diagram:  inshore 
and  offshore  profiles  in  red  and  generated  FM  T/S  profiles  in  blue  at  different  offshore  locations 
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Figure  11. 

Cross-section  plots  of  feature  model  temperature  (left)  and  salinity  (right)  off  Cape  Blanco. 
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Figure  12. 

Sampling  locations  for  the  CC  (red  'x')-  and  CUC-FMs  (black  ‘x’),  superimposed  on  the  1-deg 
Levitus  (green  'x')  temperature  at  the  surface  layer.  The  six  upwelling  sections  are  shown  by  blue 
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Figure  13.  Evolution  of  40  days  of  surface  temperature  and  velocity.  Sec  Section  6.1  for  details 
on  the  evolution. 
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Figure  14.  Evolution  of  40  days  of  salinity  and  velocity  fields  at  surface. 
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Figure  15:  Evolution  of  40  days  of  eddy  kinetie  energy  at  the  surfaee 
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Figure  16.  A  near-coastal  focus  region  (36-44N,  129-121W)  showing  a  month-long  evolution  of 
the  mesoscale  eddy  field  on  the  surface  starting  from  the  meander-upwelling  FORMS-based 
initialization. 
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Figure  17.  Evolution  of  40  days  of  temperature  and  veloeity  fields  at  300  m.  See  Seetion  6.2  for 
details  on  the  evolution. 
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Figure  18.  A  near-coastal  focus  region  (36-44N,  129-121W)  showing  a  month-long  evolution  of 
the  mesoscale  eddy  field  at  300m  starting  from  the  meander-upwelling  FORMS-based 
initialization. 
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Figure  19.  Evolution  of  40  days  of  temperature  and  velocity  fields  at  1000m.  Sec  Section  5.3  for 
details  on  the  evolution. 
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Figure  20.  A  SeaWifs  image  from  September  21,  2001  (left)  shows  eddy  activity  reasonably 
similar  to  the  model  simulation  on  day  20.  Note  the  resemblance  in  eddies  at  42N,  40N  and  36N. 
Furthermore,  the  complicated  multiple  eddies  between  37  and  40N  in  the  model  is  somewhat 
similar  to  the  smaller  eddies  visible  between  38  and  40N  in  the  image  on  the  left.  (Image 
Courtesy  of  NASA/Goddard  Space  Flight  Center,  the  SeaWiFS  Project  and  GeoEye,  Scientific 
Visualization  Studio). 
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